Electric Circuit Theory and the 
Operational Calculus 

By JOHN R. CARSON 

Note: This is the first of three installments by Mr. Carson which will 
embody material given by him in a course of lectures at the Moore School 
of Electrical Engineering, University of Pennsylvania, May, 1925. No 
effort has been spared by the author to make his treatment clear and as 
simple as the subject matter will permit. The method of presentation is 
distinctively pedagogic. To electrical engineers and to engineering in- 
structors, this exposition of the fundamentals of electric circuit theory and 
the operational calculus should be of great value. — Editor. 

Foreword 

THE following pages embody, substantially as delivered, a course 
of fifteen lectures given during the Spring of 1925 at the Moore 
School of Electrical Engineering of the University of Pennsylvania. 

After a brief introduction to the subject of electric circuit theory, 
the first chapters are devoted to a systematic and fairly complete 
exposition and critique of the Heaviside Operational Calculus, a 
remarkably direct and powerful method for the solution of the differ- 
ential equations of electric circuit theory. 

The name of Oliver Heaviside is known to engineers the world over: 
his operational calculus, however, is known to, and employed by, 
only a relatively few specialists, and this notwithstanding its remark- 
able properties and wide applicability not only to electric circuit 
theory but also to the differential equations of mathematical physics. 
In the writer's opinion this neglect is due less to the intrinsic diffi- 
culties of the subject than to unfortunate obscurities in Heaviside's 
own exposition. In the present work the operational calculus is 
made to depend on an integral equation from which the Heaviside 
Rules and Formulas are simply but rigorously deducible. It is the 
hope of the writer that this mode of approach and exposition will be 
of service in securing a wider use of the operational calculus by en- 
gineers and physicists, and a fuller and more just appreciation of 
its unique advantages. 

The second part of the present work deals with advanced problems of 
electric circuit theory, and in particular with the theory of the propaga- 
tion of current and voltage in electrical transmission systems. It is 
hoped that this part will be of interest to electrical engineers gener- 
ally because, while only a few of the results are original with the 
present work, most of the transmission theory dealt with is to be 
found only in scattered memoirs, and there accompanied by formid- 
able mathematical difficulties. 
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While the method of solution employed in the second part is largely 
that of the operational calculus, I have not hesitated to employ 
developments and extensions not to be found in Heaviside. Fcr 
example, the formulation of the problem as a Poisson integral equation 
is an original development which has proved quite useful in the actual 
numerical solution of complicated problems. The same may be said 
of the Chapter on Variable Electric Circuit Theory. 

In view of its two-fold aspect this work may therefore be regarded 
either as an exposition and development of the operational calculus 
with applications to electric circuit theory, or as a contribution to 
advanced electric circuit theory, depending on whether the reader's 
viewpoint is that of the mathematician or the engineer. 

I have not attempted in the text to give adequate reference to the 
literature of the subject, now fairly extensive. In an appendix, 
however, there is furnished a list of original papers and memoirs, for 
which, however, no claim to completeness is made. 

CHAPTER I 

The Fundamentals of Electric Circuit Theory 

While a knowledge, on the reader's part, of the elements of electric 
circuit theory will be assumed, it seems well to start with a brief 
review of the fundamental physical principles of circuit theory, the 
mode of formulating the equations, and some general theorems which 
will prove useful subsequently. 

First, the circuit elements are resistances, inductances, and con- 
densers. The network is a connected system of circuits or branches 
each of which may include resistance, inductance and capacitance 
elements together with mutual inductance, and mutual branches. 

The equations of circuit theory may be established in a number 
of different ways. For example, they may be based on Maxwell's 
dynamical theory. In accordance with this method, the network 
forms a dynamic system in which the currents play the role of veloci- 
ties. If we therefore set up the expressions for the kinetic energy, 
potential energy and dissipation, the network equations are deducible 
from general dynamic equations. 

The simplest, and for our purposes, a quite satisfactory basis for 
the equations of circuit theory are found in Kirchhoff' s Laws. These 
laws state that 

1. The total impressed force taken around any closed loop or 
circuit in the network is equal to the potential drop due to (a) resist- 
ance, (b) inductive reaction and (c) capacitive reactance. 
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2. The sum of the currents entering any branch point in the net- 
work is always zero. 

Let us now apply these laws to an elementary circuit in order to 
deduce the physical significance of the circuit elements. 

Consider an elementary circuit consisting of a resistance element 
R, an inductance element L and a capacity element C in series, and 
let an electromotive force E be applied to this circuit. If / denote the 
current in the circuit, the resistance drop is RT, the inductance drop 
is Ldl/dt and the drop across the condenser is Q/C where Q is the 
charge on the condenser. It is evident that Q and I are related by the 
equation I = dQ/dt or Q= \ldt. Now apply Kirchhoff's law relating 
to the drop around the circuit : it gives the equation 

RI+LdI/dt+Q/C=E. 

Multiply both sides by /: we get 

The right hand side is clearly the rate at which the impressed force is 
delivering energy to the circuit, while the left hand side is the rate 
at which energy is being absorbed by the circuit. The first term 
RF is the rate at which electrical energy is being converted into heat. 
Hence the resistance element may be defined as a device for con- 
verting electrical energy into heat. The second term — — LP is the 

III _ 

rate of increase of the magnetic energy. Hence the inductance 
element is a device for storing energy in the magnetic field. The 

third term — Q 2 /2C is the rate of increase of the electric, energy. 

Hence the condenser is a device for storing energy in the electric field. 

In the foregoing we have isolated and idealized the circuit elements. 
Actually, of course, every circuit element dissipates some energy in 
the form of heat and stores some energy in the magnetic field and 
some in the electric field. The analysis of the actual circuit element, 
however, into three ideal components is quite convenient and useful, 
and should lead to no misconception if properly interpreted. 

Now consider the general form of network possessing n independent 
meshes or circuits. Let us number these from 1 to n, and let the 
corresponding mesh currents be denoted by I lt J 2 • • • • In. Let 
electromotive forces E h E 2 . . . . E„ be applied to the n meshes or 
circuits respectively. Let Ljj, Rjj, Cjj denote the total inductance, 






(1) 
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resistance and capacity in series in mesh j and let Ljh, Rjk, Cjk denote 
the corresponding mutual elements between circuit j and k. Now 
write down Kirchhoff's equation for any circuit or mesh, say mesh 1; 
it is 

.... + ( Lu d i+ R u ,+ ±-fd l )l„=E 1 

Corresponding equations hold for each and every one of the n 
meshes of the network. Writing them all down, we have the system 
of equations 

{ L »i +R »+kf dt ) h+ - ■ ■ + ( L 4 l +R " + M d, ) h=E ' 
{ L 4 t +R - + hf dt ) h+ - ■ ■ + ( L «m +R -+Lf d ') I - E - 

The system of simultaneous differential equations (1) constitute 
the canonical equations of electric circuit theory. The interpreta- 
tion and solution of these equations constitute the subject of Electric 
Circuit Theory, and it is in connection with their solution that we 
find the most direct and logical introduction to the Operational Cal- 
culus. 

As an example of the appropriate mode of setting up the circuit 
equations, consider the two mesh network shown in sketch 1. Writ- 
ing down Kirschhoff's Law for meshes 1 and 2, respectively, we have 

( L ^+ R >+kf dt ) h+M P*= B > 
+ M P> + ( L 4 +R > + kf d, ) I °-= E > . 

In this case the self and mutual coefficients are given by 
Lm=L\ L 2 2 = Li Li2=Z, 2 i=+M 

Rll=Rl Ro2 = R-l 7?12 =i?21 = o 

The conventions adopted for the positive directions of currents and 
voltages are indicated by the arrows. The sign of the mutual in- 
ductance M will depend on the relative mode of winding of the two 
coils. 
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Now write down Kirchhoff's Law, or the circuital equation for 
the network of sketch 2. They are 

{ (h+u>ji +(*.+*.)+ (£+£)/*}/, 

-( L '$i +R > + kf dt ) 11 

Comparison with equations (1) shows that 

J-#ii =Li + Z,3 Z/22 = i2 + J-'3 Ll2 = Ln= — Z,3 

i?U=i?l + i?3 i?22 = -K2 + -^3 -^12=^21= — i?3 

j_ = J_ = _JL 

Cu C21 C3 



C\ 1 Ci C3 



i-i+i. 

C22 Ci Cz 



It should be observed that the signs of the mutual coefficients Rn, 
L12, C12 are a matter of convention. For example if the conventional 
directions of 1% and £2 are reversed, the signs of the mutual coefficients 
are reversed. 
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Sketch 2 

The system of equations (1) possesses two important properties 
which are largely responsible for the relative simplicity of classical 
electric circuit theory. First, the equations are linear in both currents 
and applied electromotive forces. Secondly, the coefficients Ljk, 
Rjk, Cjk are constants. Important electrotechnical problems exist, 
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in which these properties no longer obtain. The solution, however, 
for the restricted system of linear equations with constant coefficients 
is fundamental and its solution can be extended to important prob- 
lems involving non-linear relations and variable coefficients. These 
extensions will be taken up briefly in a later chapter. 

Another important property is the reciprocal relation among the 
coefficients; that is Ljk = Lkj: Rjk = Rkj, and Qk=Ckj. It is easily 
shown that these reciprocal relations mean that there are no con- 
cealed sources or sinks of energy. Again important cases exist where 
the reciprocal relations do not hold. Such exceptions, however, 
while of physical interest do not affect the mathematical methods 
of solution, to which the reciprocal relation is not essential. 

Returning to equation (1) we shall now derive the equation of 
activity. Multiply the first equation by Ii, the second by I 2 , etc. and 
add : we get 

£22t^+£2:2tj5«»+ (2) 

The right hand side is the rate at which the applied forces are supplying 
energy to the network. The first term on the left is the rate of in- 
crease of the magnetic energy 

while the second term is the rate of increase of the electric energy 

The last term, 22 Rjk Ij Ik ' is the rate at which electromagnetic 
energy is being converted into heat in the network. Consequently 
in the electrical network, the magnetic energy is a homogeneous 
quadratic function of the currents, the electric energy is a homogene- 
ous quadratic function of the charges, and the rate of dissipation 
is a homogeneous quadratic function of the currents. In Maxwell's 
dynamical theory of electrical networks, these relations were written 
down at the start and the circuit equations then derived by an ap- 
plication of Lagrange's dynamic equations to the homogeneous quad- 
ratic functions. 

Returning to equations (1), we observe that, due to the presence of 
the integral sign, they are integro-differential equations. They are, 
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however, at once reducible to differential equations by the substitu- 
tion I = dQ/dt, whence they become 

(3) 

(l^+«4+&.)e,+ • ■ +(L„4+R„4+S m )Q„ = E„. 

Here, as a matter of convenience, we have written l/Cjk—Sjk. It 
is often more convenient, at least at the outset, to deal with equations 
(3) rather than (1). 

The Exponential Solution 

In taking up the mathematical solution of equations (1), we shall 
start with the exponential solution. This is of fundamental import- 
ance, both theoretically and practically. It serves as the most 
direct introduction to the Heaviside Operational Calculus, and in 
addition furnishes the basis of the steady-state solution, or the theory 
of alternating currents. 

To derive this solution we set Ei = Fie* 1 and put all the other 
forces £2. • • E„ equal to zero. This latter restriction is a mere matter 
of convenience, and, in virtue of the linear character of the equations, 
involves no loss of generality. 

Now, corresponding to E\ = Ftf™, let us assume a solution of the 
form 

Ij = J je U (; = l,2..n) 

where Jj is a constant. So far this is a pure assumption, and its cor- 
rectness must be verified by substitution in the differential equations. 
Now if Ij = Jje* 1 , it follows at once that 

and 



fl j dt=\l j =\j/ t . 



Now substitute these relations in equations (1) and cancel the com- 
mon factor e^K We then get the system of simultaneous equations 

(\Lu+.Rii+l/XCii)/i+. . + (\L ln +R ln + l/\C ln )J n = F u 

(XL»i+.R*i+l/XC«)/i+. . + (XL 2 „+i? 2 „-r-l/XC 2 „)/„ = o J 

(4) 

OiLnl+Bnl+l/\C n i)J l +. . + (\L„ H +R HH +l/\C nn )Jn = 0. 

JNe note that this is a system of simultaneous algebraic equations 
from which the time factor has disappeared. It is this that makes 
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the exponential solution so simple, since we can immediately pass 
from differential equations to algebraic equations. In these algebraic 
equations, n in number, there are n unknown quantities J\, ■ ■ Jn-' 
These can therefore all be uniquely determined. We thus see that 
the assumed form of solution is possible. 

The notation of equations (4) may be profitably simplified as fol- 
lows : write 

\L jk + Rjk + 1 AC* = -/a (X) = z jk 



and we have 



z-nJi+ZiiJo + . .+z 2n J n = o, 



(5) 



Zn lJl + Z n 2J2 + • .-\-ZnnJn=0. 

The solution of this system of equations is 

Mji(\) Ma 

Jj= ~DW 1 = ~D ' 
and r Mjx n u Fx \t 

where D is the determinant of the coefficients, 

Zn Z12 Sl3 Zm 

Z2l Znn 223 2 2» 

Z31 Z32 Z 3" 



Znl Z n i Z nn 



(6) 



(7) 



and Mji is the cofactor, or minor with proper sign, of the jth column 
and first row. 

I shall not attempt to discuss the theory of determinants on which 
this solution is based. 1 We may note, however, one important 
property. Since Zj k = z k j, Mjk = Mkj. From this the Reciprocal 
Theorem follows immediately. This may be stated as follows: 

If a force Fe u is applied in the jth mesh, or branch, of the net- 
work, the current in the &th mesh, or branch, is by the foregoing 



I) 



Fe 



Now apply the same force in the &th mesh, or branch, then the cur- 
rent in the jth mesh is 

-pFe ■ 

1 For a remarkably concise and complete discussion of the exponential solution 
by aid of the theory of determinants, see Cisoidal Oscillations, Trans. A. I. K. E., 
1911, by G. A. Campbell. 
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Comparing these expressions and remembering that Mkj = Mjk, it 
follows that the current in the feth branch corresponding to an expo- 
nential impressed e.m.f. in the jth branch, is equal to the current in 
the jth branch corresponding to the same e.m.f. in the £th branch. 
This relation is of the greatest technical importance. 

In many important technical problems we are interested only in 
two accessible branches, such as the sending and receiving. In such 
cases, where we are not concerned with the currents in the other 
meshes or branches, it is often convenient to eliminate them from 
the equation. Thus suppose that we have electromotive forces E\ 
and Ei in meshes 1 and 2 and are concerned only with the currents 
in these meshes. If we solve equations 3, 4, . . n, « — 2 in number, 
for I3 ■ • In in terms of I\ and I* and then substitute in (1) and (2) 
we get 

ZnIi-\-Z\2li = E\, 

Znli+Znli-Et. 

The Steady State Solutions 

The steady state solution, on which the whole theory of alternating 
currents depends, is immediately derivable from the exponential 
solution. Let us suppose that E»=E 3 = . . . =E„=o and that E\ = 
F cos {oit — 6). Now by virtue of the well known formula in the 
theory of the complex variable, cos x = \e' Xj r\e~' x , we can write 

E 1 = $Fe i ( ul - B ) -\-\Fe- i ^ t - B \ 

= \ (cos d - i sin 0) Ft"*+l (ccs + * sin 0) Fe ~ iat > (9) 

= ±F'e it > t +\F"e- iut . 

Now, by virtue of this formula, the applied electromotive force E\ 
consists of two exponential forces, one varying as e 1ut and the other 
as e~ mt . Hence it is easy to see that the currents are made up of 
two components, thus 

Ij = J/c^' + J/'c-i" 1 (j = 1 ,2 . . n) (10) 

and we have merely to use the exponential solution given above, 
substituting for X.z'co and —iu respectively. That is, 

1 F' F" 

Ji' =4 ^TTR and //'H 



or 



2 Zji(ita) ' ~ Zjii-iu) 

1 Fe ~ ie M +! jjg? e -^ 

2 Zji(iu) 2 Zji(-ico) 
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The second term is the conjugate imaginary of the first, so that 

= R F c i{u,i ~ e) 
Zji(i<a) 

D F i(at-0-<t>) 

p 

cos (u)t—d—4>). 



| Z(iu>) 

We thus arrive at the rule for the steady state solution : 

If the applied e.m.f. is F cos (co/-0), substitute iw for d/dt in the 
differential equations, determine the impedance function 

Z(to)=D(to)/I(to) (11) 

by the solution of the algebraic equations, and write it in the form 

ZW= \Z(i»)\J+. (12) 

Then the required solution is 

I = -JL- cos (<o/-0-0). (13) 

I Z(t<a) | 

This in compact form contains the whole theory of the symbolic solu- 
tion of alternating current problems. 

The Complementary Solution 

So far in the solutions which we have discussed the currents are of 
the same type as the impressed forces: that is to say in physical 
language, the currents are "forced" currents and vary with time in 
precisely the same manner as do the electromotive forces. Such 
currents are, however, in general only part of the total currents. In 
addition to the forced currents we have also the characteristic oscilla- 
tions; or, in mathematical language, the complete solution must 
include both particular and complementary solutions. This may be 

shown as follows: Let iY In Devolutions of the complementary 

equations, 



{ L 4 +R » + ^J dt ) h ' + ■ ■ + { L 4i +R '" + -kf di ) L '=°' 
( L 4 +R " + kf dl )'' + ■ ■ + { L 4& R ™ + -Lf di ) I "'=°- 



(14) 



CIRCUIT THEORY AND OPERATIONAL CALCULUS 695 

Then if I\ . . . I n is a solution of (1), I\-\-I\, . . . J»-f-J»', is also a 
solution. 

To derive the solution of the complementary system of equations 
(14), assume that a solution exists of the form 

// = ///' = 1|2..») 

so that d/dt = \ and j dt = \/\. Substitute in equations (14) and 
cancel out the common factor e xt . Then we have 

Zn(A)//+ +Z ln (\)J„' = o, 

(15) 

Z»i(X)/i'+ +Z nn {\)J n ' = 0. 

This is a system of n homogeneous equations in the unknown quan- 
tities J\ ', . . Jn . The condition that a finite solution shall exist is 
that, in accordance with a well known principle of the theory of 
equations, the determinant of the coefficients shall vanish. That is, 

Zu (X) Z in (X) 



D(\) = 



Z»i(X) Znn (X) 



= 0. (16) 



Consequently the possible values of X must be such that this equation 
is satisfied. In other words, X must be a root of the equation D(X) =o. 
Let these roots be denoted by Xi, X 2 . . X m . Then, assigning to X any 
one of these values, we can determine the ratio J//Jk from any (n — 1) 
of the equations. That is to say, if we take 

Ii'-CiP+cJ*+ .... +C m e Xmt , (17) 

substitution in any (» — 1) of the equations determines iY> • • In. 
The m constants &, . . C m are so far, however, entirely arbitrary, 
and are at our disposal to satisfy imposed boundary conditions. 

This introduces us to the idea of boundary conditions which is of the 
greatest importance in circuit theory. In physical language the 
boundary conditions denote the state of the system when the electro- 
motive force is applied or when any change in the circuit constants 
occurs. The number of independent boundary conditions which 
can, in general, be satisfied is equal to the number of roots of the 
equation D(\)=o. Evidently, therefore, it is physically impossible 
to impose more boundary conditions than this. On the other hand, 
if this number of boundary conditions is not specified, the complete 
solution is indeterminate : That is to say, the problem is not correctly 
set. As an example of boundary conditions, we may specify that the 
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electromotive force is applied at time t = o, and that at this time all 
the currents in the inductances and all the charges on the condensers 
are zero. 

So far we have been following the classical theory of linear differ- 
ential equations. We have seen that the forced exponential solution 
and the derived steady state solution are extremely simple and are 
mere matters of elementary algebra. The practical difficulties in the 
classical method of solutions begin with the determination of the 
constants d, . . Cm of the complementary solution as well as the 
roots Xi, . . Xm of the equation D{\)=o. It is at this point that 
Heaviside broke with classical methods, and by considering special 
boundary conditions of great physical importance, and particular 
types of impressed forces, laid the foundations of original and powerful 
methods of solution. We shall therefore at this point follow Heavi- 
side's example and attack the problem from a different standpoint. 
In doing this we shall not at once take up an exposition of Heaviside's 
own method of attack. We shall first establish some fundamental 
theorems which are extremely powerful and will serve us as a guide 
in interpreting and rationalizing the Heaviside Operational Calculus. 

CHAPTER II 

The Solution when an Arbitrary Force is Applied to the 
Network in a State of Equilibrium 

In engineering applications of electric circuit theory there are 
three outstanding problems: 

(1) The steady state distribution of currents and potentials when 
the network is energized by a sinusoidal electromotive force. This 
problem is the subject of the theory of alternating currents which 
forms the basis of our calculations of power lines and the more elabor- 
ate networks of communication systems. 

(2) The distribution of currents and potentials in the network in 
response to an arbitrary electromotive force applied to the network 
in a state of equilibrium, i.e., applied when the currents and charges 
in the network are identically zero. 

(3) The effect on the distribution of currents and potentials of 
suddenly changing a circuit constant or connection, such as opening 
or closing a switch, while the system is energized. 

We shall base our further analysis of circuit theory on the solutions 
of problem (2), for the following reasons: 

(A) It is essentially a generalization of the Heaviside problem and 
its solution will furnish us a key to the correct understanding and 
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interpretation of operational methods and lead to an auxiliary formula 
from which the rules of the Operational Calculus are directly de- 
ducible. 

(B) The solution of problem (2) carries with it the solution of 
problem (3) and also serves as a basis for the theory of alternating 
currents. 

(C) The solution of problem (2) leads directly to an extension of 
circuit theory to the case where the network contains variable ele- 
ments: i.e., circuit elements which vary with time and in which non- 
linear relations obtain. 

Problem (2) is therefore the fundamental problem of circuit theory 
and the formula which we shall now derive may be termed the funda- 
mental formula of circuit theory. 

Consider a network in any branch of which, say branch 1, a unit 
e.m.f. is inserted at time t = o, the network having been previously in 
equilibrium. By unit e.m.f. is meant an electromotive force which 
has the value unity for all positive values of time (t^o). Let the 
resultant current in any branch, say branch n, be denoted by A, n (t). 
A n \ (/) will be termed the indicia! admittance of branch n with respect 
to branch 1 — or, more fully, the transfer indicial admittance. 

The indicial admittance, aside from its direct physical significance, 
plays a fundamental role in the mathematical theory of electric cir- 
cuits. In words, it may be defined as follows: The indicial admittance, 
A„\{t), is equal to the ratio of the current in branch n, expressed as a 
time function, to the magnitude of the steady e.m.f. suddenly inserted 
at time t = o in branch 1. It is evidently a function which is zero for 
negative values of time and approaches either zero or a steady value 
(the d.c. admittance) for all actual dissipative systems, as t approaches 
infinity. It may be noted that, aside from its mathematical determi- 
nation, which will engage our attention later, it is an experimentally 
determinable function. 

We note, in passing, an important property of the indicial admit- 
tance Ajk(t), which is deducible from the reciprocal theorem: 2 this 
is that Ajk(t) =Akj{t). That is to say, the value of the transfer 
indicial admittance is unchanged by an interchange of the driving 
point and receiving point. It is therefore immaterial in the expression 
Ajk(t) whether the e.m.f. is inserted in branch j and the current 
measured in branch k, or vice-versa. In general, unless we are con- 
cerned with particular branches, the subscripts will be omitted and 
we shall simply write A(t), it being understood that any two branches 

2 Exceptions to this relation exist where the network contains sources of energy 
such as amplifiers. These need not engage our attention here. 
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or a single branch (for the case of equal subscripts) may be under 
consideration. 

From the linear character of the network, it is evident that if a 
steady e.m.f. E = E T is inserted at time *=f,jthe network being in 
equilibrium, the resultant current is 

E T -A{t-r). 

Generalizing still further, suppose that steady e.m.fs. E , Ei, E 2 , . . . E„ 
are impressed in the same branch at the respective times t , n, t 2 
... . t«; the resultant current is evidently 

n 

E A{t)+EiA{t-T,)+ . . +EnA{t-r n ) =^EjA(t-rj). (18) 



j = o 



To apply the foregoing to our problem we suppose that there is 
applied to the network, initially in a state of equilibrium, an e.m.f. 
E{t) which has the following properties. 

1. It is identically zero for t<o. 

2. It has the value E{o) for o^t<At. 

3. It has the value E(o)+ A ,£ for A/</<2A/. 

4. It has the value E^+AxE+A^ for 2A*</<3A/. 

In other words it has the increment AjE at time t=jM. 
Evidently then the resultant current /(/) is 

EoA(t)+AiEA(t-M)+ . . +A n E.A(t-nAt). 

Now evidently if the interval At is made shorter and shorter, then 
in the limit At-*-dt and jAt = r and 

AjE = ^-E{r)dr. 

Passing to the limit in the usual manner this summation becomes a 
definite integral and we get 

I{t)=E(o)A{t)+ £A{t-r)j- T E{r)dr. (19) 

Finally by obvious transformations of the expression we arrive at the 
fundamental formula of circuit theory 

I(t) = ^-f t A(t-r)E(r)dr, (20) 

a/Jo 

= % f l E{t-T)A{r)dr. (20-a) 
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For completeness we write down the following equivalents of (20) 
and (20-a) 

1(0 -4 (*)£(*)+ j*A'(t-r)E{j)dr % (20-b) 

= A(o)E(t)+ rA'( T )E(t-T)dT, (20-c) 

= E{o)A{t)+ f'E'(t-T)A(T)dT, (20-d) 

Jo 

= E(o)A(t)+ rE'(r)A(t-T)dT. (20-e) 

where the primes denote differentiation with respect to the argu- 
ment. Thus A'{t)=d/dt A(t). 

These equations are the fundamental formulas which mathematic- 
ally relate the current to the type of applied electromotive force and 
the constants and connections of the system, and constitute the first 
part of the solution of our problem. The most important immediate 
deductions from these formulas are expressed in the following theorems. 

1. The indicial admittance of an electrical network completely 
determines, within a single integration, the behavior of the network 
to all types of applied electromotive forces. As a corollary, a knowledge 
of the indicial admittance is the sole information necessary to com- 
pletely predict the performance and characteristics of the system, 
including the steady state. 

2. The applied e.m.f. and the inidical admittance are similarly 
and coequally related to the resultant current in the network. As a 
corollary the form of the current may be modified either by changing 
the constants and connections of the network or by modifying the 
form of the applied e.m.f. 

3. Since the applied e.m.f. may be discontinuous these formulas 
determine not only the building up of the current in response to an 
applied e.m.f. but also its subsidence to equilibrium when the e.m.f. 
is removed and the network left to itself. In brief, formulas (20) 
reduce the whole problem to a determination of the indicial admittance 
of the network. In addition, as we shall see, they lead directly to 
an integral equation which determines this function. 

It is of interest to show the relation between formulas (20) and the 
usual steady state equations. To do this let the e.m.f., applied at 
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time t = o, be E sin (a>/+0). Substitution in formula (20-b) and 
rearrangement gives 

I(t)=A(o)E sin (ut+0) 

+E sin (ut+d) [ cos utA'{j)&t 
•/o 

-E cos (ut+6) f sin arA'OAdr (21) 

Jo 

where A'(t)=jA{t). 



Now this can be resolved into two parts 

E sin (u/+0) | ^ (o) +J cos urA'tfd 

-E cos (w/+0) j y sin wTi4'(T)dr f 
which is the final steady state, and 

— E sin (wt+9) I cos arA'^dr 

sin o) T A'(T)dT 



) 

T I 



(22) 



(23) 



which is the transient distortion, which ultimately dies away for suf- 
ficiently large values of time. 

To correlate the foregoing expressions for the steady state with 
the usual formulas we observe that if the symbolic impedance of the 
network at frequency oj/27t be denoted by Z(ua), and if we write 



Z(iu>) 

then the steady state current is 

E[a(w). sin (&j/-|-0)+/3(a>) . cos («*+&)]. 

Comparison with (22) gives at once 

a(u) =A(o)+ f cos wt A'Mdr, (24) 

p(to)= - (""smut A'Mdr. ■ (25) 

^0 
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The Integral Equation for the Indicia! Admittance 

So far we have tacitly assumed that the indicial admittance is 
known. As a matter of fact its determination constitutes the essential 
part of our problem. It is, in fact, the Heaviside problem, and its 
investigation, to which we now proceed, will lead us directly to the 
Operational Calculus. 

Heaviside's method in investigating this problem was intuitive and 
"experimental". We, however, shall establish a very general integral 
equation from which we shall directly deduce his methods and ex- 
tensions thereof. 

Let us suppose that an e.m.f. e"', where p is either positive real 
quantity or complex with real part positive, is suddenly impressed 
on the network at time t = o. It follows from the foregoing theory 
that the resultant current I(t) will be made up of two parts, (1) a 
forced exponential part which varies with time as e pl , and (2) a com- 
plementary part which we shall denote by y(t). The exponential or 
"forced" component is simply e p '/Z(p), where Z{p) is functionally of 
the same form as the usual symbolic or complex impedance Z(ico). 
It is gotten from the differential equations of the problem, as explained 
in a preceding section, by replacing d"/dt" by p", cancelling out the 
common factor e pl , and solving the resulting algebraic equation. The 
complementary or characteristic component, denoted by y(t), depends 
on the constants and connections of the network, and on the value of 
p. It does not, however, contain the factor e pl and it dies away for 
sufficiently large value of /, in all actual dissipative systems. Thus 

Now return to formula (20-a) and replace E(t) by e . We get 

7(/)=4-^' f'A(T)e-> T dT 
at Jo 

which can be written as 

It '! Ct " F ' A ( T ) e ~ lTdr - ct "f X ' A{r)e- pT UT } . 
Carrying out the indicated differentiation this becomes 

I(t)=pei"f A(T)e-*> T dT-pt<" f A{T)e-* T dT+A(t). (27) 
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Equating the two expressions (26) and (27) for /(/) and dividing 
through by e pl we get 

1 +y(t)e-> t =p f M A{T)e-" T dr-p f A{r)e-» T dT+A{t)e-* 1 . (28) 
Z{p) «/o Jt 

This equation is valid for all values of /. Consequently if we set 
/= so ( and if the real part of p is positive, only the first term on the 
right and the left hand side of the equation remain, the rest vanishes, 
and we get 

This is an integral equation 3 valid for all positive real values of p, 
which completely determines the indicial admittance A(t). It is on this 
equation that we shall base our discussion of operational methods and 
from which we shall derive the rules of the Operational Calculus. 
Equations (20) and (29) constitute a complete mathematical formula- 
tion of our problem, and from them the complete solution is obtainable 
without further recourse to the differential equations, or further con- 
sideration of boundary conditions. 

To summarize the preceding: we have reduced the determination 
of the current in a network in response to an electromotive force 
E(t), impressed on the network at reference time t = o, to the mathe- 
matical solution of two equations: first the integral equation 

mp)=r me - p ' dt m 

and second, the definite integral 

m =^J l A(t-r)E(r)dr. (20) 

It will be observed that in deducing these equations we have merely 

postulated (1) the linear and invariable character of the network and 

(2) the existence of an exponential solution of the type e pt /Z(p) for 

positive values of p. Consequently, while we have so far discussed 

these formulas in terms of the determination of the current in a finite 

network, they are not limited in their application to this specific 

problem. In this connection it may be well to call attention explicitly 

to the following points. 

3 An integral equation is one in which the unknown function appears under the 
sign of integration. (29) is an integral equation of the Laplace type. If Z(p) is 
specified, A(t) is uniquely determined. Methods for solving the integral equations 
are considered in detail later, in connection with the exposition of the Operational 
Calculus. The phrase "all positive values of p" will be understood as meaning all 
values of p in the right hand half of the complex plane. 
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The formulas and methods deduced above apply not <3rrly to finite 
networks, involving a finite system of linear equations, but to infinite 
networks and to transmission lines, involving infinite systems of equa- 
tions, and partial differential equations: in fact to all electrical and 
dynamical systems in which the connections and constants are linear 
and invariable. 

Secondly the variable determined by formula (20) and (29) need 
not, of course, be the current. It may equally well be the charge, 
potential drop, or any of the variables with which we may happen 
to be concerned. This fact may be explicitly recognized by writing 
the formulas as: 

im = f me - ptdt ' m 

x{t)=j t £h{t-r)E{r)dr. (31) 

Here E(t) is the applied e.m.f., x(t) is the variable which we desire to 
determine (charge, current, potential drop, etc.), and 

x=E/II(p) (32) 

is the operational equation. H(p) therefore corresponds to and is 
determined in precisely the same way as the impedance Z(p), but it 
may not have the physical significance or the dimensions of an im- 
pedance. Similarly in character and function, h(t) corresponds to the 
indicial admittance, though it may not have the same physical sig- 
nificance. It is a generalization of the indicial admittance and may be 
appropriately termed the Heaviside Function. Similarly H(p) may 
be termed the generalized impedance junction. 

CHAPTER III 

The Heaviside Problem and the Operational Equation 

The physical problem which Heaviside attacked and which led to 
his Operational Calculus was the determination of the response of a 
network or electrical system to a "unit e.m.f." (zero before, unity after 
time t = o) with, of course, the understanding that the system is in 
equilibrium when the electromotive force is applied. His problem 
is therefore, essentially that of the determination of the indicial 
admittance. In our exposition and critique of Heaviside's method of 
dealing with this problem we shall accompany an account of his own 
method of solution with a parallel solution from the corresponding 
integral equation of the problem. 
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Heaviside's first step in attacking this problem was to start with the 
differential equations, and replace the differential operator d/dt by 
the symbol p, and the operation j dt by 1/p, thus reducing the equa- 
tions to an algebraic form. He then wrote the impressed e.m.f. as 
1 (unity), thus limiting the validity of the equations to values of t = o. 
The formal solution of the algebraic equations is straightforward and 
will be written as 

h = l/H(p) (33) 

where h is the "generalized indicial admittance," or Heaviside func- 
tion (denoting current, charge, potential or any variable with which 
we are concerned) and H(p) is the corresponding generalized im- 
pedance. Thus, if we are concerned with the current in any part 
of the network, we write 

A=\/Z{p). (34) 

The more general notation is desirable, however, as indicating the 
wider applicability of the equation. 
The equations 

h = l/H(p) 

A=l/Z(p) 

are the Heaviside Operational Equations. They are, as yet, purely 
symbolic and we have still the problem of determining their explicit 
meaning and in particular the significance of the operator p. 

Comparison of the Heaviside Operational Equations with the 
integral equations (29) and (30) of the preceding chapter leads to 
the following fundamental theorem. 

The Heaviside Operational Equations 

A = 1/Z(p) 

h = l/H(P) 

are merely the symbolic or short-hand equivalents of the corresponding 
integral equations 

The integral equations, therefore, supply us with the meaning and sig- 
nificance of the operational equations, and from them the rules of the Oper- 
ational Calculus are deducible. 
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By virtue of this theorem, we have the advantage, at the outset, 
of a key to the meaning of Heaviside's operational equations, and a 
means of checking and deducing his rules of solution. This will 
serve us as a guide throughout our further study. 

Returning now to Heaviside's own point of view and method of 
attack, his reasoning may be described somewhat as follows : — 
The operational equation 

h = l/H(p) 

is the full equivalent of the differential equations of the problem and 
must therefore contain the information necessary to the solution 
provided we can determine the significance of the symbolic operator 
p. The only way of doing this, when starting with the operational 
equation, is one of induction : that is, we must compare the operational 
equation with known solutions of specific problems and thus attempt 
to infer by induction general rules for interpreting the operational 
equation and converting it into the required explicit solution. 

The Power Series Solution 

Let us start with the simplest possible problem : the current in 
response to a "unit e.m.f." in a circuit consisting of an inductance L 
in series with a resistance R. 

The differential equation of the problem is 

at 

where A is the indicial admittance. Consequently replacing d/dt by 
p, the operational equation is 

A = ~^— 
PL+K 

The explicit solution is easily derived : it is 

A=±(l-e-">) 

where a = R/L. Note that this makes the current initially zero, so 
that the equilibrium boundary condition at t = o is satisfied. 

Now suppose that we expand the operational equation in inverse 
powers of p: we get, formally, 

A=± I =±2 L_ = J_ ¥SL-.(!L\*jl.(SL\*-.(!L\4+ I 

pLl+a/p Rpi+a/p Rip K p /T V V ) 

by the Binomial Theorem. 
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Now expand the explicit solution as a power series in / : it is 

1 j at (at)* (gQ» J 

i? I 1! 2! " r 3! J ' 

Comparing the two expansions we see at once that the operational 
expansion is converted into the explicit solution by assigning to the 
symbol 1/p" the value t"/n\. It was from this kind of inductive 
inference that Heaviside arrived at his power series solution. 

Now there are several important features in the foregoing which 
require comment. In the first place the operational equation is 
converted into the explicit solution only by a particular kind of ex- 
pansion, namely an expansion in inverse powers of the operator p. 
For example, if in the operational equation 

A= l_ a/p 



R 1+a/p 

we replace 1/p by t/ll we get 

. 1 at 
A = 



R l + at 

which is incorrect. Furthermore, if we expand in ascending instead 
of descending powers of p, namely 

A=±{l-(p/a)+(p/a)*- [ 

no correlation with the explicit solution is possible and no significance 
can be attached to the expansion. We thus infer the general principle, 
and we shall find this inference to be correct, that the operational 
equation is convertible into the explicit solution only by the proper 
choice of expansion of the impedance function, or rather its reciprocal. 

In the second place we notice that in writing down the operational 
equation and then converting it into the explicit solution no con- 
sideration has been given to the question of boundary conditions. 
This is one of the great advantages of the operational method: the 
boundary conditions, provided they are those of equilibrium, are auto- 
matically taken care of. This will be illustrated in the next example : 
- Let a "unit e.m.f." be impressed on a circuit consisting of resistance 
R, inductance L, and capacity C: required the resultant charge on the 
condenser. 

The differential equation for the charge Q is 



{ i S+4 +i/c ) Qmi - '="• 
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Consequently the operational formula is 
1 



Q = 



Lp?+Rp+1/C 

1 1 U R A U 1 

where o = vand b=- 



Lpn+a/p+b/p - L —~» LC 

This can be expanded by the Binomial Theorem as 

«-i5pi»-(y+?)+(7+?) , -(7^) ,+ "f- 

Performing the indicated operations and collecting in inverse powers 
of p, the first few terms of the expansion are: — 

_L J i_£i_fi_i_f!._i_£i_fi_£i_i_ j 
Lp 2 I p p 2 ^~ p 3 ~*~ p* p b p 6 " 1- " ' ) 

where C\ = a 

c 2 = b-a 2 

c 3 = 2ab — a 3 

c 4 = 6 2 -3o 2 6-r-a 4 

c 6 = Zab 2 -4a 3 b-\-a 5 

c 6 = & 3 -6a 2 6 2 +5a 4 6-a 6 



We infer therefore that in accordance with the rule of replacing 
1/p* by t n /n\ the solution is: — 

1 j t- t 3 t* P t 6 \ 

Q = L I 2!" Cl 3!" C2 4! +C3 5! +C4 6! " ' ' J ' 

Owing to the complicated character of the coefficients in the expan- 
sion, the series cannot be recognized and summed by inspection. If, 
however, we put R = o) then a = o, and the series becomes 

whence 

Q = C{l-cos(t/VLC)\. 

We have still to verify this solution by comparison with the explicit 
solution of the differential equation. This is of the form 

Q = C+k 1 e Xil +k 2 e Xlt 



708 BELL SYSTEM TECHNICAL JOURNAL 

where k\ and k 2 are constants which must be chosen to satisfy the 
boundary conditions and X 1( X 2 are the roots of the equation 

LX 2 +i?X + l/C = o. 

Now since we have two arbitrary constants we satisfy the equilibrium 
condition by making Q and dQ/dt zero at l = o, whence 

C+ki+k 2 = o, 
\iki-}-\ 2 k 2 =o, 



and 

We have also 



fei=X 2 C/(X 1 -X 2 ), 

k 2 =\iC/(h 2 -\l). 



Writing down the power series expansion of 

<2 = C+fe ie X, '+fc 2 e X2 ', 

then 

Q = (C+k l +k 2 ) + (k 1 \ 1 +k 2 \ 2 )^ 

+ (W+foX a 2 )^ + .... 

Introducing the values of k u k 2 , \i, X 2 given above and comparing 
with the power series derived from the operational solution we see 
that they are identical term by term. 

This example illustrates two facts. First the power series expansions 
may be complicated, laborious to derive and of such form that they 
cannot be recognized and summed by inspection. In fact in arbitrary 
networks of a large number of meshes or degrees of freedom the 
evaluation of the coefficients of the power series expansion is extremely 

laborious. 

On the other hand, in such cases, the solution by the classical 
method presents difficulties far more formidable — in fact insuperable 
difficulties from a practical standpoint. First there is the location 
of the roots of the function H(\), which in arbitrary networks is a 
practical impossibility without a prohibitive amount of labor. Sec- 
ondly there is the determination of the integration constants to satisfy 
the imposed boundary conditions : a process, which, while theoretically 
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straightforward, is actually in practice extremely laborious and com- 
plicated. We note these points in passing; a more complete estimate 
of the value of the power series solution will be made later. 

To summarize the preceding: Heaviside, generalizing from specific 
examples otherwise solvable, arrived at the following rule: — 

Expand the right hand side of the operational equation 

h = l/H(p) 
in inverse powers of p : thus 

hc*a +a 1 /p+a 2 /p 2 + . . . +a n /p n + 

and then replace —j, by t n /n\. The operational equation is thereby con- 
P 

verted into the explicit power series solution : — 

h=a +a 1 t/l\+a i P/2\+ . . -\-a„t n /n\+ . . . (35) 

As stated above, this rule was arrived at by pure induction and 
generalization from the known solution of specific problems. It can- 
not, therefore, theoretically be regarded as satisfactorily established. 
The rule can, however, be directly deduced from the integral equation 



W) -jf*W^* 



pH{p) 

To its derivation from this equation we shall now proceed. 

First suppose we assume that h{t) admits of the power series ex- 
pansion 

h -\-h l t/ll-\-h 2 f i /2l + 

Substitute this assumed expansion in the integral, and integrate 
term by term. The right hand side of the integral equation becomes 
formally 

ho/p+h 1 /p* + h i /p 3 + 

by virtue of the formula 



X 



00 1" .. 1 



pt =p^+i i ^p>°- 



n\ 

Now expand the left hand side of the integral equation asymptotically 
in inverse process of p : it becomes 

a /p+a i /p i +a- i /p 3 + 

where 

a +a l /p+a 2 /p*+ .... 
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is the asymptotic expansion of I /Hip). Comparing the two ex- 
pansions and making a term by term identification, we see that 

h n = a n and 

h(t)=a +a i t/V.+a 2 t 2 /2\+ 

which agrees with the Heaviside formula. 

This procedure, however, while giving the correct result has serious 
defects from a mathematical point of view. For example, the asym- 
totic expansion of l/H(p) has usually only a limited region of con- 
vergence, and it is only in this region that term by term integration 
is legitimate. Furthermore we have assumed the possibility of ex- 
panding h(t) in a power series: an assumption to which there are 
serious theoretical objections, and which, furthermore, is not always 
justified. A more satisfactory derivation, and one which establishes 
the condition for the existence of a power series expansion, proceeds 
as follows : — 

Let 1/H(P) be a function which admits of the formal asymptotic 

expansion 

00 

^a n /P" 

o 

and let it include no component which is asymptotically representable 
by a series all of whose terms are zero, that is a function <j>(p) such 
that the limit, as p—><x>, of p n (f>(p) is zero for every value of n. Such 
a function is e~ p . With this restriction understood, start with the 
integral equation, and integrate by parts: we get 



l--=h(p)+ f" e~W»{t)dt 

(p) Jo 



H(P) 

where h (n \t) denotes d"/dt"h(t). Now let p approach infinity: in the 
limit the integral vanishes and by virtue of the asymptotic expansion 

00 

l/H(p)«>^?a n /p", (36) 

o 

1/H(p) approaches the limit a . Consequently 

h(o) =a . 
Now integrate again by parts: we get 

p(l/H(p)-a )=hW(o)+£°°e-W 2 Kt)dt. 
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Again let p approach infinity: in the limit the left hand side of the 
equation becomes fli and we have 

h^(o)= ai . 

Proceeding by successive partial integrations we thus establish the 
general relation 

h^{o)=a„. 

But by Taylor's theorem, the power series expansion of //(/) is simply 

h(t)=h(o)+h("(o)t/l\-\-hW(o)t 2 /2\+ 

whence, assuming the convergence of this expansion, we get 

h(t) =a +a 1 t/H+a 2 l n -/2l+ . . . = ^a„t"/nl (35) 

o 

which establishes the power series solution. It should be carefully 
noted, however, that it does not establish the convergence of the 
power series solution. As a matter of fact, however, I know of no 
physical problem in which II (p) satisfies the conditions for an asymp- 
totic expansion, where the power series solution is not convergent. 
On the other hand many physical problems exist, including those 
relating to transmission lines, where a power series solution is not 
derivable and does not exist. 

The process of expanding the operational equation in such a form 
as to permit of its being converted into the explicit solution is what 
Heaviside calls "algebrizing" the equation. In the case of the power 
series solution the process of algebrizing consists in expanding the 
reciprocal of the impedance function in an asymptotic series, thus 

\/II{p)«>a +a 1 , p+a i /p 2 + , 



Regarded as an expansion in the variable p, instead of as a purely 
symbolic expansion, this series has usually only a limited region of 
convergence. This fact need not bother us, however, as the series 
we are really concerned with is 

fl +a 1 //l!+a 2 / 2 /2!+ 

It is interesting to note in passing that the latter series is what Borel, 
the French mathematician, calls the associated function of the former, 
and is extensively employed by him in his researches on the summa- 
bility of divergent series. 

The process of "algebrizing," as in the examples discussed above, 
may often be effected by a straight forward binomial expansion. 
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In other cases the form of the generalized impedance function H(p) 
will indicate by inspection the appropriate procedure. A general 
process, applicable in all cases where a power series exists, is as follows. 
Write 

l/H(p) = l/H(±-)=G(x). (36) 

Now expand G(x) as a Taylor's series : thus formally 

GW=G( ) + CW(o)^+G( 2 '( )^+ . . . 
where 

Denote - , by a„, replace x n by l/p u , and we have 
G(x) = l/H(p)=a +a 1 /p+a 2 /p 2 + . . . 

This process of "algebrizing" is formally straightforward and always 
possible. As implied above, however, in many problems much shorter 
modes of expansion suggest themselves from the form of the function 

H(p). 

We note here, in passing, that the necessary and sufficient condi- 
tions for the existence of a power series solution is the possibility of 
the formal expansion of G(x) as a power series in x. 

At this point a brief critical estimate of the scope and value of the 
power series solution may be in order. As stated above, in a certain 
important class of problems relating to transmission lines, a power 
series does not exist, though a closely related series in fractional 
powers of / may often be derived. Consequently the power series 
solution is of restricted applicability. Where, however, a power 
series does exist, in directness and simplicity of derivation it is superior 
to any other form of solution. Its chief defect, and a very serious 
defect indeed, is that except where the power series can be recog- 
nized and summed, it is usually practically useless for computation 
and interpretation except for relatively small values of the time /. 
This disadvantage is inherent and attaches to all power series solu- 
tions. For this reason I think Heaviside overestimated the value 
of power series as practical or working solutions, and that some of 
his strictures against orthodox mathematicians and their solutions 
may be justly urged against the power series solution. He was quite 
right in insisting that a solution must be capable of either interpre- 
tation or computation and quite right in ridiculing those formal 
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solutions which actually conceal rather than reveal the significance 
of the original differential equations of the problem. On the other 
hand, the following remark of his indicates to me that Heaviside 
has a quite exaggerated idea of the value and fundamental character 
of power series in general: "I regret that the result should be so 
complicated. But the only alternatives are other equivalent infinite 
series, or else a definite integral which is of no use until it is evalu- 
ated, when the result must be the series (135), or an equivalent one." 
As a matter of fact the properties of most of the important functions 
of mathematical physics have been investigated and their values 
computed by methods other than series expansions. I may add that 
in technical work the power series solution has proved to be of re- 
stricted utility, while definite integrals, which Heaviside * particu- 
larly despised, have proved quite useful. 

The Expansion Theorem Solution 5 

We pass now to the consideration of another extremely important 
form of solution. Heaviside gives this solution without proof: we 
shall therefore merely state'the solution and then derive it from the 
integral equation. 

Given the operational equation 

h = l/H(p) 

which has the significance discussed above: i.e., the response of the 
network to a "unit e.m.f.". The explicit solution may be written as 

h= m+Xfrmptj (38) 

where p u pi . . . p n are the n roots of the equation 

II(p)=o 
and 

*<*-[£?&],-* (39) 

As remarked above, this solution, referred to by him as The Ex- 
pansion Theorem, was stated by Heaviside without proof; how he 
arrived at it will probably always remain a matter of conjecture. 
Its derivation from the integral equation is, however, a relatively 
simple matter, though in special cases troublesome questions arise. 

4 Vide a remark of his to the effect that some mathematicians took refuge in a 
definite integral and called that a solution. 

6 This terminology is due to Heaviside. A more appropriate and physically 
significant expression would be "The Solution in terms of normal or characteristic 
vibrations." 
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The derivation of the expansion solution from the integral equation 

follows immediately from the partial fraction expansion 

_J_ = _L_ ■ V J (40) 

pll(p) pH{oy j^{p-pj)PjH'(pj) 

where pi, pi . . . pn are the roots of the equation II(p) =o, and 

//'fe)=||/w} f=? , («) 

Partial fraction expansions of this type are fully discussed in treatises 
on algebra and the calculus and the conditions for their existence 
established. Before discussing the restrictions imposed on H(p) by 
this expansion, we shall first, assuming its existence, derive the ex- 
pansion theorem solution. 

By virtue of (40) the integral equation is 

The expansion on the left hand side suggests a corresponding expan- 
sion on the right hand side: that is, we suppose that 

h(t)=h (t)+h 1 (t)+fa(t)+ . . .+hn(t) (43) 

and specify that these component functions shall satisfy the equations 

pm'X'w** (44) 

(p-ttow f™*** j=u '-- M - (45) 

It follows at once from (43) and direct addition of equations (44) 
and (45) that (42) is satisfied and hence is solved provided h , . . h n 
can be evaluated from (44) and (45). 
Now since 

f" e»e-*tdt—r^r (46) 

Jo p — h 
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provided the real part of X is not positive (a condition satisfied in all 
network problems), we see at once that equations (42) and (43) are 
satisfied by taking 

Ut) = ho = j^ y (47) 

Consequently from (43) and (47) it follows that 

which establishes the Expansion Theorem Solution. 

As implied above, the partial fraction expansion (40), on which 
the expansion theorem solution depends, imposes certain restrictions 
on the impedance function H(p). Among these are that H(p) must 
have no zero root, no repeated roots, and \/II{p) must be a proper 
fraction. In all finite networks these conditions are satisfied, or by 
a slight modification, the operational equation can be reduced to 
the required form. The case of repeated roots, which may occur 
where the network involves a unilateral source of energy such as an 
amplifier, can be dealt with by assuming unequal roots and then 
letting the roots approach equality as a limit. Without entering 
upon these questions in detail, however, we can very simply and 
directly establish the proposition that the expansion theorem gives 
the solution whenever a solution in terms of normal or characteristic 
vibrations exists. The proof of this proposition proceeds as follows. 

It is known from the elementary theory of linear differential equa- 
tions that the general solution of the set of differential equations, 
of which the operational equation is h = l/H(p), is of the form 



h{t)=Co+^C j e Pit 



where pj is the jth root of H(p) =o, and C, C\ . . . C„ are constants of 
integration which must be so chosen as to satisfy the system of dif- 
ferential equations and the imposed boundary conditions. The 
summation is extended over all the roots of H(p) which is supposed 
not to have a zero root or repeated roots. 
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Now substitute this known form of solution in the integral equation 
of the problem and carry out the integration term by term. We get 

Setting p = o, we have at once 

Co = l/H(o). (50) 

To determine Q let p = pj+q where q is a small quantity ultimately 
to be set equal to zero, and write the equation as 

If now p=pj-\-q and q approaches zero, this becomes in the limit 

pjH'(pj)Cj = l (52) 

or 

c -» (53) 

whence W-^g+^gWK ™ 

which is the Expansion Theorem Solution. 

We shall not attempt to discuss here cases where the expansion 
solution breaks down though such cases exist. In every such case, 
however, the breakdown is due to the failure of the impedance func- 
tion H(p) to satisfy the conditions necessary for the partial fraction 
expansion (40), and correlatively the non-existence of a solution in 
normal vibrations. Furthermore, it is usually possible by simple 
modification to deduce a modified expansion solution. It may be 
added here, that while the proof given above is also limited implic- 
itly to finite networks, the expansion solution is valid in most trans- 
mission line problems. 

Let us now illustrate how the expansion solution works by applying 
it to a few simple examples. Take first the case considered in the 
preceding chapter in connection with the power series solution. Re- 
quired the charge Q on a condenser C in series with an inductance L 
and resistance R in response to a "unit e.m.f." The operational 
equation is 

Q= Lp 2 +Rp+l/C 
n 1 1 

° r V L p 2 + 2ap+a> 2 
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where a = R/2L and u 2 = l/LC. 

The roots of the equation H(p) =o are the roots of the equation 

p 2 + 2ap+u>- = o 
whence 

pi = — a-\-\/a 2 — o> 2 = — a-\-fi, 

p2= —a— v a 2 — a) 2 = — a — /3. 

Also H'{p) =2L(p+a), so that 

H'(p x )=2QL 

H'(p 2 ) = -20L 
and 

l/H(o) = l/Lu 2 = C. 

Inserting these expressions in the Expansion Theorem Solution 
(38), we get 

n = r e —( effl e ' '\ 
V 2/3LW-/3 a+p)' 

It is now easy to verify the fact that this solution satisfies the differ- 
ential equations and the boundary condition Q = o and dQ/dt = o at 
time t = o. 

If cc>a, /3 is a pure imaginary 

= iay/l-(a/u) 2 = Ua' 
and 

e~ al a/ cos a>'/ + Qf sin co'/ 



Q = C- 



u'L a 2 -\-co' 2 



In connection with this problem we note two advantages of the 
expansion solution, as compared with the power series solution: (1) 
it is much simpler to derive from the operational equation, and (2) 
its numerical computation is enormously easier. A table of expo- 
nential and trigometric functions enables us to evaluate Q for any 
value of / almost at once whereas in the case of the power series solu- 
tion the labor of computation for large values of / is very great. A 
third and very important advantage of the expansion solution in this 
particular problem is that without detailed computation we can 
deduce by mere inspection the general character of the function 
and the effect of the circuit parameters on its form : an advantage 
which never attaches to the power series solution. 

This last property of the particular solution above is extremely 
important. The ideal form of solution, particularly in technical 
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problems, is one which permits us to infer the general character and 
properties of the function and the effect of the circuit constants on 
its form, without detailed solutions. A solution which possesses 
these properties, even if its exact computation is not possible without 
prohibitive labor, is far superior to a solution which, while com- 
pletely computable, tells us nothing without detailed computation. 
It is for this reason that some of the derived forms of solution, dis- 
cussed later, are of such importance. In fact a solution which re- 
quires detailed computation before it yields the information implied 
in it is merely equivalent to an experimentally determined solution. 
Unfortunately the advantages attaching to the expansion solution 
of the specific problem just discussed, do not, in general, characterize 
the expansion solution. The following disadvantages should be 
noted. First, the location of the roots of the impedance function 
H(p) is practically impossible in the case of arbitrary networks of 
more than a few degrees of freedom. In the second place, when the 
number of degrees of freedom is large it is not only impossible to 
deduce the significance of the solution by inspection, but the com- 
putation becomes extremely laborious. In such cases, the practical 
value of the expansion solution depends, just as in the power series 
solution, on the possibility of recognizing and summing the expan- 
sion. This will be clear in the case of transmission lines, where the 
roots of H(p) are infinite in number and the direct computation of 
the expansion solution (except in the case of the non-inductive cable) 
is quite impossible. 

CHAPTER IV 

Some General Formulas and Theorems for the Solution 
of Operational Equations 

We have seen that the operational equation 

h = l/H(p) 
is the symbolic or short-hand equivalent of the integral equation 

and from the latter we have deduced two very important forms of the 
Heaviside solution. In recognizing the equivalence of these two 
equations we have a very great advantage and are able, in fact, to 
base the Operational Calculus on deductive instead of inductive 
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reasoning. In this chapter we shall employ this equivalence to estab- 
lish certain general formulas and theorems for the solution of oper- 
ational equations. That is to say, we shall make use of the principles 
that (1) any method applicable to the solution of the integral equation 
supplies us with a corresponding method for the solution of the 
operational equation, and (2) a solution of any specific integral equa- 
tion gives at once the solution of the corresponding operational 
equation. We turn therefore to a brief discussion of the appropriate 
methods for solving the integral equation. 

It may be said at the outset, that the solution of the integral equa- 
tion, like the evaluation of integrals, is a matter of considerable art 
and experience; in other words there is not, in general, a straight- 
forward procedure corresponding to the process of differentiation. 

On the other hand, as a purely mathematical question, it is always 
possible to invert the integral equation and write down h(t) as an 
explicit function in the form of an infinite integral. For example 
it may be shown from the Fourier Integral that 



*»-!/.' 



! a(o>) . , 
sin tco.au 



where a(a>) is defined by 

T«U-a(«)+»0(«). 

U(iu) 

Later on we shall briefly consider the Fourier Integral; for the 
present the preceding formula will not be considered further. In 
certain problems it is of value; for the explicit derivation of h(t), 
however, it is usually too complicated to be of any use except in the 
hands of professional mathematicians. As a matter of fact, a direct 
attack on this formula would be equivalent to abandoning the unique 
simplicity and advantages of the whole Operational Calculus. 

It has been noted above that any solution of the integral equation 
supplies a solution of the corresponding operational equation. This 
principle enables us to take advantage of the fact that a very large 
number of infinite integrals of the type 



X 



f(t)e-"'dt 



have been evaluated. The evaluation of every infinite integral of this 

type supplies us, therefore, with the solution of an operational equation. 

Of course, not all the operational equations so solvable have physical 

significance. Many, however, do. Below is a list of infinite integrals 
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with their known solutions, accompanied by the corresponding 
operational equation and its explicit solution. All of these solutions 
are directly applicable to important technical problems. It may be 
remarked in passing that the infinite integrals have for the most 
part been evaluated by advanced mathematical methods which need 
not concern us here. 

Table of Infinite Integrals, the Corresponding Operational Equations, 
and Their Explicit Solutions 

1 



V fo 



e -pt e -\tdt= 



h=-^-=e-K 



(b) f °e-i> t - 1 dt = l/p n + l , 
«/o nl 

X 00 1 1 

V-xt Vp 

h=Vp=VVvi. 

(d) f 

h = 



e - P t. ( 2 0" dt _ 1 



o 1.3.5... (2«-l)Vr/ p»Vp' 

Vp (2/)" 1 



p» 1.3.5... (2n-l) Vtt/* 
/" 1 

n\ e (p+X)** 1 ' 



(e) / e-t"—.e- Xt dt^ 
«/o 



(0 



h = P = !%-x, 

(p-\-\) n + 1 nl ' 

r e - Pl Jh^U l=e -^u 

tVt' 



h=pe~ 2 ^P = J- 



-X/< g-2V\/> 



r°° p-w p -2s\p 
h = Vpe-^^=— 7 =. 

y/irt 
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X 



(h) / e-i> l s\n\tdt = 



,2_L\2' 



'o £ 2 +X 



(i) I e~* % cos X/ d/ = 
Jo 



£ 2 +X 2 ' 



*-^» =cosW - 



G) /"V 



.-* cos \tdt= ' ' 



(p+ri'+X" 



(k) J fr^r^d.»Hf (H .^. +v . 



h = - 7 JL==j {u). 
vV+x 2 



(m) / e-"'J (Vt 2 -^)dt= . --, 

J\ vp 2 +i 

h = -J=e~^' m =olort<\ 
Vp 2 + 1 

= / (V/ 2 -X 2 ) for/>X. 
(n) fj e - t " J ^t)dt = y( L ^y r 2 =/> 2 + X 2 , 

(p) f X e-P'e*I {\t)dt= .-L- _ -, 
Vl+2X/p 
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In formulas (1), (m), (n), J„{x) denotes the Bessel function of order 
n and argument x. In formula (p), I {x) denotes the Bessel function 
J (ix) where i=\/—l. 

This list might be greatly extended. As it is, we are in possession 
of a set of solutions of operational equations which occur in important 
technical problems and which will be employed later. 

The foregoing emphasize the practical and theoretical importance 
of recognizing the equivalence of the integral and operational equa- 
tions. With this equivalence in mind, the solution of an operational 
equation is often reduced to a mere reference to a table of infinite 
integrals. Heaviside did not recognize this equivalence. As a 
consequence many of his solutions of transmission line problems are 
extremely laborious and involved and in the end unsatisfactory 
because expressed in involved power series. 

Not all the infinite integrals corresponding to the operational 
equations of physical problems have been evaluated or can be recog- 
nized without transformation. This statement corresponds exactly 
with the fact that a table of integrals is not always sufficient but 
must be supplemented by general methods of integration. We turn, 
therefore, to stating and discussing some general Theorems applicable 
to the solution of Operational Equations. 

In the derivation of the operational theorems, which constitute the 
general rules of the Operational Calculus, the following proposi- 
tion, due to Borel and known as Borel's theorem, will be frequently 
employed.* 

// the functions f(t),fi(t), andf 2 (t) are defined by the integral equations 



F(P) = f X 'f(t)e-P'dt 

Fi(P)= rh{t)t-m 

•A) 

f*(p)= f~MQrm 



and if the functions F, F\ and F 2 satisfy the relation 

F(p) = F 1 (p).F 2 (p) 

* For a proof of this important theorem the reader is referred to Borel, " Lecons 
sur les Series Divergentes" (1901), p. 104; to Bromwich, "Theory of Infinite Series," 
pp. 280-281; or to Ford, "Studies on Divergent Series and Summability," pp. 93-94 
(being Vol. II of the Michigan University Science Series, published by Macmillan). 
The proof depends on Jacobi's transformation of a double integral: see Edward's 
"Integral Calculus," 1922, Vol. II, pp. 14-15. 
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then 



f(0= f'Mr)Mt-r)dr 
J(> 

= f'h{r)h{t-r)dr. 
Jo 



The operational theorems will now be stated and briefly proved 
from the integral equation identity. 

Theorem I 

If in the Operational Equation 

h = l/II(p) 

the generalized impedance function H(p) can be expanded in a sum of 
terms, thus 

1 1,1. 1 



II{P) Ih(p) ' H,(P) ■ '"^Hnipy 
and if the auxiliary operational equations 

1 



lh = 

h„_ = 



Hi(p) 

1 



Ih(P) 

can be solved, then 

//=//l + //2 + . . • +kn. 

This theorem is too obvious to require detailed proof: in fact it is 
self evident. The power series and expansion theorem solutions are 
examples of its application. In general, however, the appropriate 
form of expansion of 1/H(p) will depend on the particular problem 
in hand. The theorem, as it stands is a formal statement of the fact 
that solutions can often be obtained by an appropriate expansion 
whereas the equation cannot be solved as it stands. 

Theorem II 
If h = h{t) and g=g(t) are defined by the operational equations 
h = l/H(p) 
g = l/pH(p) 



then 



g(t)= f'h{r)dr. 
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To prove this theorem we start with the integral equations 

pkp) = £ h{t)e ' m ' 

The second of these is in form for an immediate application of Borel's 
theorem since 

1 _ 1^ 1 
p*H(p) 'p'pH(py 

The functions /i and / 2 of Borel's theorem then satisfy the equations 

\= rfi(t)e-*'dt, 

p «/o 

It follows at once that 

Mt)=h(t) 

whence by Borel's theorem 

g{t)=£h{r)dr. 

Theorem III 

If h=h(t) and g=g(t) are defined by the operational equations 

h = l/H(p) 
g = P/H(p) 



then 



««-5*0 



provided h{o) =o. 

The integral equations of the problem are 
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Integrating the first of these by parts we have, 

mm -T h M + jf h ' (t)e ~ m 

where h'{t)=d/dth{t). 

If h{6) =o, we have at once 

Comparison with the integral equation for g{t) shows at once that 
g{t)=h'{t), since the integral equation determines the function 
uniquely. 

Theorems II and III establish the characteristic Heaviside Opera- 
tions of replacing 1/p by / dt and p by d/dt. 

Jo 

Theorem IV 
If in the operational equation 

h = l/II(p) 
the generalized impedance function can be factored in the form 

II(p)=II 1 (p).H 2 (p) 
and if the auxiliary operational equations 
hi = l /Ih(p) 
h, = l/H 2 (p) 
define the auxiliary variables hi and h 2 , then 

h(t)=~f Q 'h 1 (r)h 2 (t-r)dr 

= it£ i,2(T)i,i(t ~ T)dT - 

This theorem is immediately deducible from Borel's theorem and 
theorems II and III, as follows. 
The integral equations are 

pikprf l '> (t)e - f ' dL 
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Now define an auxiliary function g(t) by the operational equation 
1 

g ~pii(py 

Then 

and by Borel's theorem 

g(0= f l hMh*{t-T)dT 



= f h 2 (r)hi(t-T)dT. 
Jo 



From this equation it follows that g(o) =o, and hence comparing the 
operational equations for h and g, we have by aid of Theorem III 

and hence 

Jk(0-4 Ai(r)fc 8 (*-r)<*T 

4r*i(rM-r)dr. 

This theorem is extremely important, although not stated or 
employed by Heaviside himself. We shall make use of it in estab- 
lishing two important general theorems and shall have frequent 
occasion to employ it in specific problems occurring in connection 
with the subsequent discussion of transmission theory. 

Theorem V 
If h = h(t) and g = g(t) are defined by the operational equations 

1 



h = 



H(p) 

1 



where X is a positive real parameter, then 

g(f) = (l+\£dt)e-*h(f). 
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To prove this theorem we start with the integral equations 






In the first of these equations replace the symbol p by g+X : we get 

-4t • m L = f* h{t)e-*e-«dt 
g+X H{q+\) Jo 

and then to preserve our original notation replace the symbol q by £, 
whence 

(7+xjA(?+x)=X"" (/)e " x ^"*- < a > 

The integral equation in g{t) can be written as 

/ x \ i r °° 

( i+ 7) ^+x)/^+x) = y * (w - (b » 

Comparing equations (a) and (b) it follows at once from theorems 
I and II that 

From the foregoing, the following auxiliary theorem is immedi- 
ately deducible. 



Theorem Va 

If h = h{t) and g = g(t) are defined by the operational equations 

1 



h = 



5 = 



II{p) 

P 



{p+\)H{p+\) 
then 

g(t)=h(t)e-K 

The proof of this theorem will be left as an exercise to the reader. 
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Theorem VI 
If h = h(t) and g = g(t) are defined by the operational equations 
h = l/II(P) 
g=l/II(\p) 
where X is a positive real parameter, then 
g(t)=h(t/\). 
We start with the integral equations 

and in the first of these equations we replace p by \q and / by t/X, 
whence it becomes 

Now replacing the symbols q and t by p and / respectively, we have 

whence by comparison with the integral equation in g(t) it follows 
at once that 

g(t)=h(t/\). 

This theorem is often useful in making a convenient change in the 
time scale and eliminating superfluous constants. 

Theorem VII 

If h = h(t) and g = g(t) are defined by the operational equations 
1 



h = 



H(p) 



g H{p) 
where X is a positive real quantity, then 
g(t)=oior t<\ 

= h(t-\) for />X. 
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This is a very important theorem in connection with transmission 
line problems where retardation, due to finite velocity of propagation, 
occurs. Its proof proceeds as follows: 

If the auxiliary function k = k(t) is defined by the operational 
equation 

then by Theorem IV, 

g{t)=~f\{r)h{t-r)dr. (a) 

Now, corresponding to the operational equation k = e we have 
the integral equation 

e -\P 



— = / k(t)e-*'dt. 
P «/o 



The solution of this integral equation, which is easily verified by 
direct substitution in the infinite integral, is 

k{t)=o for /<X 

= 1 for />X. 

Hence equation (a) becomes 

g(f)=o for /<X 

=4l A(*- r)dr for fe\ 

dtJ\ 

= h(t-\) for f£\. 

Theorem IV, employed in the preceding proof, as stated above, is 
extremely important and we shall have frequent occasion to employ 
it in specific problems. We shall now apply it to deduce an important 
theorem which extends the operational calculus to arbitrary impressed 
forces, whereas heretofore the operational equation h = l/H(p) applied 
only to the case of a "unit e.m.f." impressed on the system. 

It will be recalled from a previous chapter that if x(t) denotes the 
response of a network to an arbitrary force /(/), impressed at time 
t = o, and if h(t) denotes the corresponding response to a "unit e.m.f.," 
then 

x(t)=j t £h(r)f{t-r)dr (31^ 



and 



fiH(p) 



= I h{t)e-t"dt. (30) 
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Now/(/) may be of such form that the infinite integral 



£ 



/o 

can be evaluated and has the value F(p)/p: thus 

f~ f{t)e-»dt = ±-F{p). (55) 

Jo P 

This is possible, of course, for many important types of applied forces, 
including the sinusoidal. 

It follows at once from Theorem IV that x(t) satisfies and is de- 
termined by the integral equation 

iM= r<t)e^dt. (56) 

p H{p) Jo 

We have thus succeeded, by virtue of Theorem IV in expressing the 
response of a network to an arbitrary e.m.f. impressed at time t = o, 
by an integral equation of the same form as that expressing the 
response to a "unit e.m.f." That is to say we have, at least formally, 
extended the operational calculus explicitly to the case of arbitrary 
impressed forces. 

We now translate the foregoing into the corresponding Operational 
Theorem. 

Theorem VIII 

If the operational equation 

h = l/H(p) 

expresses the response of a network to a "unit e.m.f." and if an arbitrary 
e.m.f. E impressed at time t = o,is expressible by the operational equation 

E=V{p) 
or the infinite integral 

f; E ^m=m 

then the response x of the network to the arbitrary force is given by the 
operational equation 

V(p) 

mpy 

and x{t) is determined by the integral equation 
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Theorem IX 
If the operational equation 

h = l/II(p) 
is reducible to the form 

F(p) 



fc = 



l+\K(p) 



where X is a real parameter, and if the auxiliary functions /=/(/) and 
k = k(t) are defined by the auxiliary operational equations 

f = F(p) 

k=K(p) 

then h(t) is determined by the Poissan Integral equation 

h{t)=f{t)-\ f l h{ T )k{t-T)dT. 

This theorem is of considerable practical importance in connection 
with the approximate and numerical solution of operational equations 
when the operational equation and the equivalent Laplace integral 
equation prove refractory. In such cases, as will be shown later, 
the numerical solution of the Poissan integral equations can often 
be rapidly and accurately effected, and in many cases the quali- 
tative properties of h{t) can be deduced from it without detailed 
numerical solution. 

The proof of this theorem proceeds as follows: 

By virtue of the relation h = l/II(p) the operational equation 



/,= "<» 



can be written as 



l+\K(p) 



*+*f$=TO 
*=TO-*f§ 



A direct application of Borel's theorem or Theorem IV gives at once 
the explicit equivalent 



W) =Rt)~\ fh{j)k(t-T)dT. 
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The preceding theorems, together with the power series and ex- 
pansion theorem solutions formulate the most important rules of the 
operational calculus, and are constantly employed in the solution of 
the electrotechnical problems. On the other hand, the table of 
infinite integrals furnishes the solution of a set of operational equa- 
tions, which are of the greatest usefulness in the systematic study of 
propagation phenomena in transmission systems which will engage 
our attention. Before taking up this study, however, we shall first 
solve a few specific problems which will serve as an introduction to 
asymptotic and divergent solutions involving Heaviside's so-called 
"fractional differentiation." 

Problem A : Current Entering the Non-inductive Cable 

The non-inductive cable is a smooth line with distributed resistance 
R and capacity C per unit length ; for the present we neglect induct- 
ance and leakage. A consideration of cable problems leads to some 
of the most interesting questions relating to operational methods, 
particularly to questions regarding divergent expansions. It would 
seem best to allow specific problems to serve as an introduction to 
these general questions. 

The differential equations of the cable are 

RI=-~-V 

dx 

(57) 

dt dx 

where x is the distance, measured along the cable from any fixed 
point, 2" is the current at point x, and V the corresponding potential. 
Replacing d/dt by the operator p, we have 

RI=-^-V 
dx 

(58) 
PCV=-I-L 



Eliminating, successively, V and / from these equations, we get 

dx 2 

and 

_9_ 2 
dx* 



PRCI = ^-J 



P RCV=£- 2 V. 
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These equations have the general solutions 

V-V*-v+Vj* (59) 

i =4^ [vie ' yx - v ^ x] < 6o > 

where 

y = VjRC. (61) 

The term in e~ yx represents the direct wave and the term in e yx the 
reflected wave. V\ and V~2 are constants which must be so chosen 
as to satisfy the imposed boundary conditions at the terminals of 
the cable. 

For the present we shall assume that the line is infinitely long so 
that the reflected wave is absent. We shall also assume that a voltage 
E is impressed directly on the cable at x = o: we have then, 

V=Ee~ x ^pcK = Ee-^ (62) 

I = yj^Ee- x ^^R = ^Ee-^ (63) 

where a denotes x 2 RC. 

To convert these to operational equations let us suppose that E 
is a "unit e.m.f." (zero before, unity after time t = o). We have 
then, in operational notation 

I = ^e-^P. (65) 

Now suppose that x = o so that a = o, in other words consider a 
point at the cable terminals. Then 

(66) 



/ 



=# 



The first of these equations means that V is simply the impressed 
voltage, zero before, unity after time t = o, as of course, it should be 
from physical considerations. 

Corresponding to the operational equation 

(66) 



-M 
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we have the integral equation 

yli3p-f™<- m - (67) 

The solution of this is known (see formula (c) of the preceding table 
of. integrals) : it is 

I = X [?^=J^. (68) 

Heaviside arrived at this solution from considering the known 
solution of the same problem in the theory of heat flow. He there- 
fore inferred that the operational equation 

has the explicit solution 

This is correct; we, however, have derived it directly from the integral 
equation of the problem and the known integral 

i r>* (69 ) 

We then see from the foregoing that, if a "unit e.m.f." is impressed 
on the cable terminals, the current entering the cable is i nitially 
infinite and dies away in accordance with the formula y/CJ-wRt. 
The case is, of course, idealized and the infinite initial value of the 
current results from our ignoring the distributed inductance of the 
cable, which, no matter how small, keeps the initial current finite, 
as we shall see later. 

Now let us go a step farther; suppose that in addition to distributed 
resistance R and capacity C, the cable also has distributed leakage 
G per unit length. The differential equations are now 

RI= --^-V 
■dx 

a (70) 

(Cp+G)V=-^ x I. 

Consequently it follows that in the operational equation for the current 
entering the cable we need only replace Cp by Cp-\-G. Therefore, 
when leakage is included, equation (66) is to be replaced by 



/ = V R 
where \ = G/C. 



*T^-JlV*+* C7i) 
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The corresponding integral equation is, of course, 

Vx^-r^w^* (72) 

We shall give two solutions of this problem; first the solution of 
the integral equation, and second the typical Heaviside solution 
directly from the operational equation. 

Equation (72) may be written as 

Now suppose that /(/) is the solution of the equation 

X — = r* J(t)e-**dt (74) 

it follows at once from Theorems (I) and (II) of the preceding chapter 
that 

/(/) = \lf ( 1+x jT*) /w ' (75) 

Also from formula (c) of the table of integrals and Theorem (Va) the 
solution of (74) is 

P -\t 
AO— =7= (76) 

whence 

™-&& +x £$*\- (77) 

The integral appearing in (77) can not be evaluated in finite terms; 
it is easily expressible as a series, however, by repeated integration 
by parts. Thus 

1^-7= dt = 2 / e-^'dVT = 2VTe-*+2\ ler^y/Tdt. 

J -V t •'0 J 

Proceeding in this way by repeated partial integration we get for the 
integral term of (77) 

*•!>*{ i+»{ + S$ + ... I. (78) 

' 1.6 1.6.0 ) 
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The straightforward Heaviside solution is obtained by expanding 
the operational equation as follows : 

-JiM^ng) , ^g) , -]^ : 

Identifying -\/p with l/y/iri (from known solutions of allied problems) 
and substituting for l/p" multiple integrations of the rath order we get 

t r^ii_J 2X/ ) (2X/) 2 , 1.3(2X0' ) (m 

It can be verified that this solution is convergent and equivalent 
to (77). 

This problem, while simple and of minor technical interest, will serve 
to introduce us to the very important and interesting question of 
asymptotic series solutions. 

An asymptotic series, for our purposes, may be defined as a series 
expansion of a function, which, while divergent, may be used for 
numerical computation, and which exhibits the behavior of the func- 
tion for sufficiently large values of the argument. 

Let us return to equation (77). We observe that the series solu- 
tion (78) of the definite integral becomes increasingly laborious to 
compute as the value of / increases. This remark applies with even 
greater force to the Heaviside solution (79) on account of the alter- 
nating character of the series. Right here we have an excellent 
example of what I regard as Heaviside's exaggerated sense of the 
importance of series solutions as compared with definite integrals. 
Consider the solution in the form of (77) as compared with Heavi- 
side's series solution (79). The former is incomparably easier to 
interpret and to compute, either by numerical integration or by 
means of an integraph or planimeter. In fact the series (79) is prac- 
tically unmanageable except for small values of t. 

Returning to the question of an asymptotic expansion of the solu- 
tion (77), we observe that the definite integral appearing in that 
equation can be written as, 

C . VT dt= f° VT dt S VT dl (80) 



CIRCUIT THEORY AND OPERATIONAL CALCULUS 72,7 

provided X is positive, as it is in thi s case. Now the value of the 
infinite integral is known; it is \/ir/\. Consequently 

furthermore, 

•°° «-x< 1 /"° i .1 p-\t i /.■» *-x< 






y/ 1 *y/t 2X^0 ty/t 

Integrating again by parts we get 

1 e-* 1 1 g~ x< 1.3 f°°_e^_ H 
*VT 2X 2 /V //" 2W0 *vf 

Continuing this process, we get 



dt. 



s; 



e~^ e-* 1 p 1_ 1.3 1.3.5 

VT XVTL 2X/ + (2X/) 2 (2X/) 3 

+ .. +( -D-ai s #=a] 



(82) 



(2X/)» 
(-1)" 1.3.5 . . . (2w+l) Z* 00 e-* 



I _! 

•/1 <■- 



= d*. 



X 2(2X)» Ji i»+V/ 

Now this series is divergent, that is, if we continue out far enough 
in the series the terms begin to increase in value without limit. On 
the other hand, if we stop with the nth term the error is represented 
by the integral term in (82) and this is less than 

i-iy 1.3.5... (2rc-l) 

Consequently the error committed in stopping with any term in the 
series is less than the value of that term. Therefore if we stop with 
the smallest term in the series, the error is less than the smallest term 
and decreases with increasing values of t. 
We can therefore write the solution (77) as 

SR + \rrRt e (2Xf (2X/) 2+ (2X/) 3 I (84) 

The first term, since \ = G/C, is simply \/G/R, the d.c. admittance 
of the leaky cable. The divergent series shows how the current 
approaches this final steady value. 
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In this particular problem no asymptotic solution is derivable 
directly from the operational equation, at least by the straight- 
forward Heaviside processes. Asymptotic solutions, however, con- 
stitute a large and important part of Heaviside's transmission line 
solutions. We shall therefore discuss next a problem for which 
Heaviside obtained both convergent and divergent series expansions. 

Problem B: Terminal Voltage on Cable with "Unit E.M.F." Impressed 
on Cable Through Condenser 

We now take up a problem for which Heaviside obtained a divergent 
solution, and which will introduce us to the theory of his divergent 
solutions and so-called "fractional differentiation." We suppose a 
"unit e.m.f." impressed on an infinitely long cable of distributed 
resistance R and capacity C per unit length through a condenser of 
capacity C : required the voltage V at the cable terminals. The 
operational equation of the problem is derived as follows: — 

We know from the problem just discussed that the current entering 
the cable whose terminal voltage is V, is, in operational notation 



V 



R 



But the current flowing into the condenser is 

C p(l-V) 

since the voltage across the condenser is 1 — V. Equating these two 
expressions we get 

V= ^= (85) 

pC +VpC/R 

which is the operational equation of the problem. 

This may be written as 

1 



V= 



1_ \CJ_ (85) 

1 



1 + Va/p' 
where 

\ r a=-rVCjR. 



CIRCUIT THEORY AND OPERATIONAL CALCULUS 739 

Now expanding this by the binomial theorem 



la a a la fa\ 2 



-"7+(f)'+"" 

'at .(at)*, , . 

-l+j-j + -2j-+ • • • • (oo) 

/2a/ (2a/) 2 (2a/) 3 \ 1 

A.1 + 1.3. + 1.3.5. + "Jy^al 

by the usual Heaviside rules of "algebrizing." 

It is worth while verifying this from the integral equation of the 
problem. We have 

-T- ^7=P " f" V(t)e-*'dt. (87) 

P 1 + Va/p -Ai w 

The left hand side can be written as 

_1 1_ Ya_ 

p — a p — a\p 

and by the formulas and theorems given in a preceding section the 
solution can be recognized at once as: — 

V(t) -0« -J %<" /"^ dr. (88) 

> 7T Jo V r 

This can also be written as 

If the definite integral of (88) is evaluated by successive partial 
integrations it will be found in agreement with the Heaviside solution 
(86). 

Now the solution (86) is in powers of / and while absolutely con- 
vergent becomes progressively more difficult to interpret and com- 
pute as the value of / increases. From (89), however, we can derive 
a divergent or asymptotic solution applicable both for interpreta- 
tion and computation, when the value of / is sufficiently large. As 



740 BELL SYSTEM TECHNICAL JOURNAL 

in the example discussed before, the asymptotic expansion results 
from repeated partial integrations; thus 

r<» p-ar 1 /*°° 1 

e- al ___ J_ P" £^_ dr 
™ a\/l 2a ^< t\/t 



_k 3 /* 

2 W, 



+ Mf ^-dr 



and finally 



fl \/^ 2aHy/t 2W ( tV^ 

^771 2a/ + (2a/)' (2a/)* + ' ' " I' ^ 

The series (90) is divergent just as is (82) of a preceding problem 
and the error committed by stopping with the smallest term, is of 
the same character and subject to the same discussion. With this 
understanding we write the solution (89) as 

w-:?yi-£+&-s + ---i- (9i) 

For large values of t (at>5) this series is accurately and rapidly 
computable. Furthermore it shows by mere inspection the be- 
havior of 7(0 for large values of /, and that it ultimately approaches 
zero as l/s/irat. 

Let us now see how Heaviside attacked this problem and how he 
arrived at a divergent solution from the operational formula. Re- 
turning to the operational equation (85), it can be written as 

V= ^Eg . (92) 

l + Vpfa 
Now expand the denominator by the binomial theorem: we get 
formally 

-(i+(£) ! H-(ir + 4 
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Heaviside's procedure at this point was as remarkable as it was suc- 
cessful. He first discarded the second series in integral powers of p 
as meaningless. He then identified y/p with l/y/irt and replaced 
p" by d"/dt" in the first series, getting 






(94) 



\firat 
or, carrying out the indicated differentiation, 

\/™t\ 2ar(2at)* (2o0* / 

which agrees with (91). 

This is a typical example of a Heaviside divergent solution for 
which he offered no explanation and no proof other than its practical 
success. His procedure in this respect is quite unsatisfactory and in 
particular his discarding an entire series without explanation is in- 
tellectually repugnant. We shall leave these questions for the present, 
however; later we shall make a systematic study of his divergent 
solutions and rationalize them in a satisfactory manner. First, 
however, we shall take up a specific problem for which Heaviside 
obtains a divergent solution without discarding any terms. 

Problem C: Current Entering a Line of Distributed L, R and C 

Consider a transmission line of distributed inductance L, resistance 
R, and capacity C per unit length. The differential equations of 
current and voltage are 

A a (95) 

dt dx 



Replacing d/dt by p, we get 

( P L+R)I=-^V 

(96) 



Equations (96) correspond exactly with (58) for the non-inductive 
cable: except that we must replace R by pL+R. For the infinitely 
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long line, therefore, the operational formula for the current entering 
the line is 

r-J P c y (97) 

where V is the voltage at the line terminals. If this is a "unit e.m.f." 
we have, as our operational equation, 



/ = 



PC— (98) 



which can be written as 



\ pL+R 



Z-j£-=i=™ (99) 



"Mi V.1+2X/P 



where \ = R/2L. 

The corresponding integral equation is 

\f-7 / * =1 e- pt Kt)dt. (100) 

From either equation (99) or (100) and formula (p) of the table of 
integrals, we see at once that the solution is 

r=^e-*i (\t) (ioi) 

where I o (A0 is the Bessel function 7 (iX/), where i = y/— 1. (The 
function is, however, a pure real.) 

Heaviside's procedure, in the absence of any correlation between 
the operational equation and the infinite integral, was quite different. 
Remarking, with reference to equation (99), that "the suggestion to 
employ the binomial theorem is obvious," he expands it in the form 

r _j£J 1 -^+t§^M , -i^5^V+ r ( 102 ) 

t-^T r p + 2i\p) 3! \p) +••) 

and replaces 1/p" by t"/n in accordance with the rule discussed in 
preceding sections. The explicit solution is then 

/=^I 1 - w+ ^-i-W + ...| a* 

a convergent solution in rising powers of /. As yet, however, he does 
not recognize this series as the power series expansion of (101), which 
it is. He does, however, recognize the practical impossibility of 
using it for computing for large values of /, and remarks "But the 
binomial theorem furnishes another way of expanding the operator 
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(operational equation), viz. in rising powers of p." Thus, returning 
to (99), it can be written as, 

\£ Vl+/>/2X 

Now expand the denominator by the binomial theorem : we get 

\L I 4X^2! UxJ 3! U\) T " f \ 2X 

He now identifies -\/p/2\ with l/\/2irX/ and replaces />" in the series 
by d"/W, thus getting finally 

i= JJ-L= \ i+J-+iW +<IM!+ I . (ice) 

^ L v/2irX/ ' ^8X/^2!(8X/) 2 ^3!(8X/) 3 ^ ' - " S 

This series solution is divergent: Heaviside recognizes it, however, 
as the asymptotic expansion of the function e~ Xt I (\t), and thus 
arrives at the solution 



I=^~e->«I (\t) (101) 

which we have obtained from our tables of integrals. 

Now the divergent expansion (106) is the well known asymptotic 
expansion of the function e I„(\t), which is usually derived by diffi- 
cult and intricate processes. The directness and simplicity with 
which Heaviside derives it is extraordinary. 

We note in this example that no integral powers of p appear in the 
divergent expansion : consequently no terms are discarded. Other- 
wise Heaviside's process is as startling and remarkable as in the 
example discussed in the preceding section. 

We shall later encounter many problems in which asymptotic 
solutions are derivable as in the preceding example. We have suffi- 
cient data, however, in these two typical examples to take up a 
systematic discussion of the theory of Heaviside's divergent solution 
of the operational equation. 

CHAPTER V 

The Theory of the Asymptotic Solution of Operational 

Equations 

A study of Heaviside's methods, as exemplified in the preceding 
examples and in many problems dealt with in his Electromagnetic 
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Theory, Vol. II, shows that they may be divided into two classes: 
(I) those of which the operational equation is of the form 

h = F(p)Vp (I) 

and (II) those of which the operational equation is of the form 

h = 4>(p k Vp) (II) 

where k is an integer. 

Heaviside himself does not distinguish between the two classes, 
but employs the following rule for obtaining asymptotic expansion 
solutions: 

// the operational equation 

h = l/H(p) 

can be expanded in the form 

h = a +a 1 p+a 2 p 2 -\- . . +a n p n + . ■ • 

(b +b 1 p+hp! i + • • +b»p n + • ■ • )Vp, (107) 

a solution, usually divergent, is obtained by discarding the first_ expansion 
entirely, except for the leading constant terms a a , replacing y/p by X/y/vt 
and p n by d"/dt n in the second expansion, whence an explicit series 
solution results. 

»-+^b+»«*+"-)73 (108) 

■ =-+^-4+ fe (yHw + ---)- < io9 > 

It should be expressly understood that Heaviside nowhere himself 
states this rule formally. He does not distinguish between the two 
cases where integral series in p do and do not appear, although very 
important mathematical distinctions are involved. Furthermore, 
in one case he modifies his usual procedure by adding an extra term 
(Elm. Th. Vol. II, pg. 42^4). It certainly represents, however, 
his usual procedure in a very large number of problems. 

A completely satisfactory theory of the Heaviside Rule, just stated, 
has not yet been arrived at although we can always verify the diver- 
gent solutions in specific problems. Furthermore, it is not as yet 
known just how general it is, though it certainly works successfully 
in a large number of physical problems to which it has been applied. 
Finally we know nothing in general as to the asymptotic character 
of the resulting expansion. In some cases it leads. to an expansion 
in which the error is less than the last term included, in others re- 
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markably enough the expansion is everywhere convergent, while 
in yet others its application leads to a series which is meaningless 
for a certain range of values of /. 

Heaviside himself gives no information which would serve us as a 
guide in informing us when the rule is applicable and when it is not. 
Consequently it becomes a matter of practical importance, not only 
to investigate the underlying mathematical philosophy of the rule 
and to establish it on the basis of orthodox mathematics, but also to 
develop if possible a criterion of its applicability. In this investiga- 
tion we shall have recourse to the integral equation of the problem. 

We shall take up first the type of problem (Class I) in which the 
operational equation is 

h =ik) =F(pw * (110) 

and assume that F{p) admits of the formal power series expansion 

F(p)=b +b 1 p+b 2 p 2 +b 3 p*+ . . . (Ill) 

The corresponding integral equation is 

Qjfi- f" h(t)e-P'dt. (112) 

■\/p Jo 

We now assume the existence of an auxiliary function k(t), defined 
and determined by the auxiliary integral equation 

F(p)= r*k(t)e-<"dt. (113) 

Jo 

Now since 

-WV>-4= (in) 

■y/p Jo v irt 

it follows from (112), (113), and (114) and Borel's Theorem, or 
Theorem IV, that 

h(t) = -U f'-^L dr. (115) 

Now if we differentiate (113) repeatedly with respect to p and put 
p = o, it follows from the expansion (III) that 

b n = (-i)" T-MOdt. (116) 

Jo w! 

This equation presupposes, it should be noted, the convergence of 
the infinite integrals for all values of n, and therefore imposes severe 
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restrictions on k{t) and hence on F(p). We shall suppose that these 
restrictions are satisfied, and discuss them later. 
Now (115) can be written as: — 

h(t) = -±= PdrMr) (1-r//)" 1 / 2 . (117) 

It can be shown that, if k(t) satisfies the restrictions underlying 
(116), the integral (117) has an asymptotic solution obtained as 
follows :— Expand the factor (1-t//)" 1/2 by the binomial theorem, 
replace the upper limit of integration by », and integrate term by 
term : thus 



m ^ t \£ km+ U x r^ dt 






l Vfr»* + -l 



Finally from (116) we get 

•»-^{fc-^+*#-*W + -} (,19) 

which agrees exactly with the Heaviside rule for this case. 

The foregoing says nothing regarding the asymptotic character 
of the solution. It is easy to see qualitatively, however, that (118) 
•and therefore (119) does represent the behavior of the definite in- 
tegral (117) for large values of t, provided k(t) converges with suffi- 
cient rapidity. 

The foregoing analysis may now be summarized in the following 
proposition : 

If the operational equation h = l/II(p) is reducible to the form 

h = F(p)Vp 
and if F{p) admits of power series expansion in p : thus 

F(p)=b +b 1 p+b 2 p*+ . . . +b„p n + 

50 that, formally, 

h = (b +b 1 p + bip 2 + . . . +b n p w + • ■ ■ )VJ 

an explicit series solution, usually asymptotic, is obtained by replacing 
Vpby l/y/vt and p n {n integral) by d"/dt n , whence 



1 /, . 1 , , 1.3 , 1.3.5 \ 



y/ivt 
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provided the function k = k(t), defined by the operational equation k = F(p) , 
and the infinite integrals . . 

f°°t»k(t)dt (w — 1,2, ) 

Jo 

exist. 

We shall now apply the foregoing theory to a physical problem 

discussed in the last section : namely, the current entering an infinitely 

long line of inductance L, resistance R and capacity C per unit length. 

It will be recalled (see equation (100) ) that the integral equation of 

this problem is 

Jg . 1 == f™ e-»I(t)dt 
\ L y/p* + 2\p Jo 

where \ = R/2L, and that the solution is 



= yj^'I (\t). 



We can derive the solution in another form appropriate for our pur- 
poses by writing 



j C 1 1 /.- 



dt 
' *» y/py/p-\-'l\ Jo 

Now since 



dt 
y/p Jo y/tct 

and 



/to Jo 



-2\t 



L= = f 

>+2X Jo 



y/p+2\ Jo y/tct 

it follows from Borel's theorem that 

* L IT Jo y/jyjt-j 

Now subject this definite integral (omitting the factor y/C/L) to 
the same process applied to (117) : we get 

ttV t ' Jo V / M J o 1 ! 

+ (»)Vo 2! g «+••• r 



748 BELL SYSTEM TECHNICAL JOURNAL 

The infinite integrals are known and have been evaluated. Sub- 
stituting their values this series becomes : — 

_J_ r 1 1 2 .3 2 1 2 .3 2 .5 2 | 

VMUl 8X* + 2!(8X0 2 + 3!(8X/) 3 + ' ' } 

which is in fact the well known asymptotic expansion of the function 
e- Xt I (\t). 

A second example may be worth while. Consider the case of an 
e.m.f. e~ x ' impressed at time t = o on a cable of distributed resistance 
R and capacity C: required the current entering the cable. The 
required formula is fl 



yivRdtJo Vt 



by obvious transformations. 

Asymptotic expansion of the definite integral as in the preceding 
example gives the asymptotic formula 



-V 



C i 1 . 1.3 , 1.3.5 



i 



icRt \2\r (2\/) 2 ' (2X/) 3 ' " f 
The operational formula of the problem is 

Applying the Heaviside Rule, we get the asymptotic expansion 

\ irRt I 2\r (2X/) 8 ^ (2X0 3 ' " » 

which agrees with the preceding formula, derived from the definite 

integral. 

We shall now discuss a specific problem in which the Heaviside Rule 
breaks down. For example let us take the preceding problem, and 

6 The derivation of the formulas in this problem is left as an exercise for the 
reader. 
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replace the applied e.m.f. e~ xt by sin cot. The formula corresponding 
to (120) is now 

I = coJ± / -7r=dr. 121 

If we now attempt to expand the definite integral of (121) in the 
same way as that of (120), we find that the process breaks down because 
each component of the infinite integral is now itself infinite. In fact 
no asymptotic solution of this problem exists. 

Let us, however, start with the operational formula : since 






e~ pl sin cot. dt=- 



it is 

Now expand this in accordance with the Heaviside Rule: we get, 
operationally, 

and explicitly 

T-- EJ1.M5.J. I 

VirRt 1 2cat (2w*) 3 ' ' I 

which is quite incorrect. 7 The incorrectness of the result will be evident 
when we remember that the final value of the current is the steady-state 
current in response to sin cot, or 



4 



^(cos co/+sin cot). (122) 

This result can be derived directly from (121) by writing it as 
\~C I f" c os cot f 1 sin cot | 

I= " \s ( cos w v ^7T dt+sln w v Vt dt s ' (123) 

If the time is made indefinitely great the upper limits of the integrals 
may be replaced by infinity. The infinite integrals are known: sub- 
stitution of their known values gives (122). 

This example illustrates the care which must be used in applying 
Heaviside's rules for obtaining divergent solutions and the importance 

7 While this series is incorrect as an asymptotic expansion of the current it has 
important significance, as we shall see, in connection with the building up of alter- 
nating currents. 
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of having a method of checking the correctness of his processes and 
results. 

We now take up the discussion of the asymptotic expansion solu- 
tions of operational equations of the type 

h=4>{p k \f£) (k integral). (123) 

In this discussion we shall, as a matter of convenience, assume that 
k = o, so that the equation reduces to the form 

h = <f>(Vp). (123a) 

This will involve no loss of essential generality, since the analytical 
theory of the two equations is precisely the same. 

The Heaviside Rule for this type of operational equation may be 
formulated as follows : 

If the operational equation h = l/H(p) is reducible to the form 

h=<i>(p k V¥) 

and if <$> admits of power series expansion in the argument, thus 

h=a +a 1 p k Vp+a 2 p 2k + i +a 3 p 2k + i v / p + ■ . . 

a series solution, usually divergent and asymptotic, is obtained by dis- 
carding integral powers of p, and writing 

h=a +(a 1 p k -{-a 3 p Sk ^-\-a s ,p^+ 2 -{- . . . )Vp- 

The explicit series solution then results from replacing -y/p by l/y/irt, 
and p n by d"/dt", whence 

rf 3fe+l dok+2 



, { d* , d"** 1 . d°*+ a . \ 



txi a - 



V T*l 

(-!)»/ 1.3... (2fe-l) n 1.3 .. . (6H-1) , \ 



VSV 1 (2/)* " J (2ft3*+i 

The theory of this series solution will be based on the following 

proposition, deducible from the identity / / dt = \/\/p. 
If the function F(p) of the integral equation 

F(p)=JJf(t)e-P'dt 

approaches 1/ y/p as p approaches zero, then /(ft ultimately behaves as 
l/y/irf. that is, if F{p)— >l/\/p as p-^o, then f(t)^>l/-\/irt as t— »oo, 
provided that /(ft converges to zero, and contains no term or factor which 
is ultimately oscillatory. 
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To illustrate what this condition means suppose that 

f ( t\ — a -4- - COS ^ 
\Zirt y/vt 
then 



f(t)e-P>dt-+a/Vp as p-*0, 



and the oscillatory term in f(t) converges to a higher order. The 
presence of such oscillatory terms vitiate, therefore, the Heaviside 
Rule : in the following discussion we shall assume that they are absent. 
We are now prepared to discuss the operational equation 

h=<t>(p k Vp) 

and for convenience shall assume that k = so that the operational 
equation becomes 

of which the corresponding or equivalent integral equation is 

jct>(Vpj=j Q h{t)e-"'dt. (123b) 

We assume that ${\/p) admits of formal power series expansion in the 
argument: thus 

<t>(Vp)=a +ai\/p+a2p+a 3 p\/p~+a i p 2 + . . . 

without, however, implying anything regarding the convergence of 
this expansion. 

We now introduce the series of auxiliary functions, g,gi,gi,g 3 

defined by the following scheme 

2(0 =//(/) -a 

gi(t)=g(0--~^ 
Virt 

g 2 (/)=/gi(0+ 1 



\f> (12 3c) 

gM=t, g . i {t)+^^L-_ 
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Successive substitutions in the integral equation (123b) and repeated 
differentiations with respect to p, lead to the set of formulas, 



f* g(t)e->'dt ~ -^ as p-*0 

Jq ' -yp 

f" Lg^e-t'dt** -^=as p->0 
»/o 2Vp 

jT^/.^Oe-^^^-^as p->0 

J /tea 
' t.g 3 (t) 




(123d) 



-m<°-^^i»o 



Now assuming that h(t) satisfies the restrictions stated in the pre- 
ceding proposition, it follows from that proposition, that 

g(t) n ai/y/rt as /— »oo 

(123e) 



2tVirt 



fA 1.3 at 

1.3.5 a 7 



From the set equations (123d) and (123e) it follows by successive 
substitutions that 

1 / 1 . 1.3 1.3.5 v 

which agrees with the series gotten by applying the Heaviside Rule. 

The defect of this derivation, which, however, appears to be in- 
herent, is that it requires us to know or assume at the outset that h(t) 
satisfies the required restrictions. Consequently an automatic ap- 
plication of the Heaviside Rule may or may not give correct results. 
On the other hand if we know that an expansion solution in inverse 
fractional powers of / exists, the Heaviside Rule gives the series with 
extraordinary directness and simplicity. 

The type of expansion solution just discussed will now be illustrated 
by some specific problems. The first problem is that of the propagated 
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voltage in the non-inductive cable in response to a "unit e.m.f". It 
will be recalled that in a preceding chapter we derived the operational 
formula 

V=e~^ (124) 

where a = x 2 RC, for the voltage at distance x from the terminal of a non- 
inductive cable of distributed resistance R and capacity C, in response 
to a "unit e.m.f." impressed at point x = 0. Heaviside's solution of 
this operational equation proceeds as follows: 

Expansion of the exponential function in the usual power series gives 

V- 1 V"? I °& ctp\ / ~ap (ap) 2 _ 

1! " r 2! 3! "*" 4! 

which may be rearranged as 

Heaviside then discards the series in integral powers of p entirely, 
replaces s/p by l/^/irt and p" by d n /dt n in the first series, and then gets 

This solution is correct, as will be shown subsequently. 

A rather remarkable feature of this solution — a point on which 
Heaviside makes no comment — is that it is absolutely convergent. 
In other words, a process of expansion which in other problems leads 
to a divergent or asymptotic solution, here results in a convergent 
series expansion. 

To verify this solution we start with the corresponding integral 
equation of the problem 

i«"V*« f" V{t)e~ pi dt. (128) 

Jo 



or 



1 -./- 
P 

It follows from this formula and theorem (V) that 



•w-jf 



4>{t)dt 
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where <f>(t) is determined by the integral equation 

Now from formula (f) of the table of integrals 

g-tt/41 



2 \xJo /V/ 



2 \ Wo /v/ 

whence 



♦»-^f 



g-O/41 



and finally ;*, 

F(/) = -4- f' ^-7= <* T - where '' = 4/ /"- (129 > 

V7T I/O TVT 

To convert this to the form of (127) we write 

FC0-- W —^ dT -^L ~l= iT - (130) 

V 7T *^° rV s/tJ' t\/t 

The value of the infinite integral is known to be unity so that 

7=1 Lr°£^lrf T . (131) 

■S/irJi' tS/t 

Now in the integral term of (131) expand e~ 1/r in the usual expo- 
nential power series and then integrate term by term : the series solu- 
tion (127) results. This series, while absolutely convergent, is difficult 
to compute for small values of /; an asymptotic expansion, which can 
be employed for computation for small values of / is gotten as follows : — 

Write (129) as 

Repeated partial integrations of this type lead to the series 

F~J5r«'{l-(4)+L8(4) , -..f. (132) 

It is interesting to note, in passing, that an asymptotic solution of 
this type does not appear to be directly deducible from the operational 
equation. We observe also that, in this problem, the series in inverse 



CIRCUIT THEORY AND OPERATIONAL CALCULUS 755 

powers of t is convergent while the series in ascending powers of / is 
divergent : the converse is the case in the problems discussed previously. 
A second specific problem may be stated as follows: 

Let a "unit e.m.f." be impressed on an infinitely long non-inductive 
cable of distributed resistance R and capacity C per unit length through 
a terminal resistance R : required the voltage V on the cable terminals. 
The formulation of the operational equation of this problem is very 
simple. It will be recalled that the operational formula for the current 
entering the cable with terminal voltage V is V\/Cp/R. But the 
current is clearly also equal to (1 — V)/R : equating these expressions 
we get 

1-V 



whence 



j<o =VVpC/R 



V=—F= — (133) 

V>A+1 



where \/y/\ = R y/ C/R. This is the required operational formula. 

To derive the Heaviside divergent expansion, expand (133) by the 
binomial theorem : thus 

v=x-Vp/\ +{p/\)-(p/\y 2 + ... 

= l-(l+£/X + (/>/X) 2 + . . WpJ\ (134) 

+ (PA + (P/A) 2 +(P/X) 3 + ...)• 

Discard the second series in integral powers of p; replace s/p by 1/x/irt 
and p n by d"/dt n in the first series, thus getting 



^ i -( i 44+^l+-)vfe 



= i- 



. *rX/ 
1 /. 1 1.3 



(135) 



iO-R+afaji---) < 136 > 



VttXA 2X/ (2X0 s 

which is the asymptotic solution of the problem. 

To verify this solution we shall consider the more general opera- 
tional equation 

1 

h= P n Vp+i ( Minte s ral ) ( 137 ) 

a form of equation to which a number of fairly important problems 
is reducible. (The parameter X of equation (133) can be eliminated 
from explicit consideration by means of theorem VI.) 
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Multiplying numerator and denominator of equation (137) by 
p H y/p—l, it becomes 

*-^^= ? fcv^-^I (138) 

and by direct partial fraction expansion, this is equivalent to 

2n+lZ* p-p m 2n+l^p-p 



»-^S£™2,-fe <m 



m=o 

where 



pm = e *»+i (m = 0,1,2. . . 2n). 
Write, for convenience, 



2n 



and consider the operational equation 

hn = _l_ (J#l v^_ Jul.) . (wo) 

" m 2n+l\p-p„V p-pj 

By the rules of the operational calculus, fully discussed in preceding 
chapters, the solution of this is 

We have now to distinguish two cases: (1) when the real part of p m 
is positive, and (2) when the real part is negative. 
Taking up case (1) first, the preceding can be written 

/i_i_ pm ' i&l f°° e-*"» T | 

I I ' \/irJo \/t T f 

= 2^+l( B** e ** C~r^!. dT ) ' ' ^ 

Repeated integration by parts of the definite integral leads to an 
asymptotic series, identical with that obtained by applying the Heavi- 
side Rule to the operational equation (137). 
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If, on the other hand, the real part of p m is negative, we write (141) as 

I I l — e Pmt \ 

*-W-5i+i| + p^_ p e pmT dr . (145) 

The term e t>m ' ultimately dies away, and the definite integral can be 
expanded asymptotically in accordance with the theory discussed 
under Rule I, again leading to an asymptotic series identical with that 
given by direct application of the Heaviside Rule to the operational 
equation. 

Consequently since the operational equation in h„ can be asymptotic- 
ally expanded by means of the Heaviside Rule, the operational equa- 
tion in h='Sh, n is similarly asymptotically expandible, and the 
Heaviside Rule is verified for equation (133). 

We have now covered, more or less completely, the theoretical rules 
and principles of the operational calculus in so far as they can be 
formulated in general terms. We shall now apply these principles 
and rules to the solution of important technical problems relating to 
the propagation of current and voltage along lines. In doing, so, while 
we shall take advantage of our table of integrals with the corresponding 
solutions of the operational equation, we shall also sketch Heavi- 
side's own methods of solution. 

We shall close this discussion of divergent and asymptotic expan- 
sions with a general expansion solution of considerable theoretical 
and practical importance in the problem of the building-up of alter- 
nating currents. It will be recalled from Theorem III that the response 
of a network of generalized operational impedance H(p) to an e.m.f. 
E{t) impressed at time t — o is given by the operational formula 

V(P) 

H{p) 

where E= V(p) is the operational equation of the applied e.m.f.: that 
is, analytically 



irV(p)= rE{t)e-»dt. 



Now suppose that the impressed e.m.f. is sin ut: then by formula (h) 
of the table of integrals 

V(P) =^ 2 (146) 
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and denoting x by x s 

*•-?+* m- (l47) 

If, on the other hand, the impressed e.m.f. is cos wt, then by formula (i) 
and 



^)-sb (148) 



£ 2 I 

Now let us consider the operational expansion suggested by the 
Heaviside processes : 

X,m "a\ l+ J) mF) 

and 

Now let us identify \/H(p) with h{t) and replace p n by d"/dt" : we get 
j l d l d 3 . I d 6 f , , A 

and 

J l d 2 l d 4 . I d 6 I . , rt , 1KQ , 

We have now to inquire into the significance of equations (152) and 
(153), derived from the operational equations of the response of the 
system of an e.m.f. sin wt and cos wt respectively, impressed at time 
/ = 0. From the mode of derivation of these expansions from the 
operational equations it might be inferred that they are the divergent 
of asymptotic expansions of the operational equations (147) and 
(149). This would certainly not be an unreasonable inference in the 
light of the Heaviside expansions we have just been considering. This 
inference is however, not correct: on the other hand, the series (152) 
and (153) have a definite physical significance, as we shall now show 
from the explicit equations of the problem. 
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By equation (31), the explicit equation for x s , given operationally 
by (147), is 

Xs= dtJ sin l0T - h (t-T)dT=J sin a)(t-T)h'(T)dT+h(o) sin at (154) 
where h'(/) =d/dt h(t). By a well known trigonometric formula, this is 

.r s = sin wt I cos ut.h'(t)dt— cos cot I sin ut.h'(t)dt-\-h(o)sm ut. 
Writing 

/»/ /too /»oo 

dt= dt- dt 

t'lis becomes 

/•oo «oo 

*.s = sin ut I cos o>t.h'(t)dt — cos «/ / sin ut.h'(t)dt 
Jo Jo 

/oo 
sin u>(t-T)h'(T)dT. (155) 

The first three terms are simply the steady-state response to the 
impressed e.m.f. sin u>/: that is, they represent the ultimate steady 
state value of x s when the transient oscillations have died away. The 
last term, which we shall denote by T s , represents the transient oscilla- 
tions which are set up when the e.m.f. is applied. Thus 

/oo 
sin <ji{t- T )h'{T)dT. (156) 

Now from (156) 

T s =-- f" V(rU.cos~a>(r-0 
and integrating by parts 

Ts = llt h{t)+ lf t Cos W (r-0;gA(r)<*T. (157) 

Repeating the process of partial integration, we get: 

1 J I n 00 73 

Ti ^H m "^l ^^-^ 3 h(r)dr. (158) 

Repeating the process again 
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This process can be repeated indefinitely, and we get 

r-( l d iii+lii- (-D"- 1 ^- 1 ^^ 

Is ~Wdt co 3 dt 3 ^co 5 dP "^ co 2 "- 1 df*"' 1 ) W 

+-^l ™^-»dWi h{T)dT - (159) 

The series expansion (159), except for the remainder term, is identical 
with the series expansion (152) derived directly from the operational 
equation. This series may be either convergent or divergent, de- 
pending on the frequency co/2ir and the character of the indicial ad- 
mittance function h(t). In the important problems of the building-up 
of alternating currents in cables and lines we shall see that, even when 
divergent, the series is of an asymptotic character and can be employed 
for computation. 

We thus arrive at the following theorem : 

If an e.m.f . sin cot is impressed at time / = on a network or system 
of generalized indicial admittance h(t), and if the transient distortion, 
T s , is defined as the instantaneous difference between the actual re- 
sponse of the system and the steady-state response, then T s can be 
expressed as the series 

( id i d* l # _ (-D- 1 &rl\ m am 
WTri?d73+»*dp ••• + «*•-! dP*-i) HW {im 

with a remainder term 



,2" 



* /'°o d 2 " +1 

J sin w(t-/) ^ 2b+1 h{r)dT. 



If the impressed e.m.f. is cos cot, the corresponding series for the trans- 
ient distortion, T c , is 

(I d? 1 d, 1 d« (-V nd2n \h(t) eun\ 

{^dJ-^di+^di*- ■■ SSTdi*)*® ( 161 ) 

with a remainder term 

(-D"f , A d*»+\ 



1) M /*°° d in+1 

{jpj cos co(r-/)^+I^(r)rfr. 



The second part of this theorem, relating to the transient distortion, 
T c , in response to an e.m.f. cos cot, is derived from formula (31) by 
processes precisely analagous to those employed above in deriving 
the series expansion for T s . The derivation will be left to the reader. 

To summarize the preceding discussion of the divergent solution of 
operational equations, it may be said that the theory is as yet rather 
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unsatisfactory. To the physicist it is unsatisfactory because he 
requires an automatic rule giving a correct asymptotic expansion by 
purely algebraic operations without investigations of remainder terms 
or auxiliary functions. Furthermore, the precise sense in which the 
expansion asymptotically represents the solution cannot be stated in 
general, but requires an independent investigation in the case of each 
individual problem. 

On the other hand when an asymptotic expansion is known to exist, 
the Heaviside Rule finds this expansion with incomparable directness 
and simplicity, the problem of justifying the expansion being a purely 
mathematical one, which usually need not trouble the physicist. 
Furthermore, on the purely mathematical side, the Heaviside Rule 
is of large interest and should lead to interesting developments in the 
theory of asymptotic expansions. 

{To be continued) 



